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Abstract. An overview of some recent results on the geometry of partial 
differential equations in application to integrable systems is given. Lagrangian 
and Hamiltonian formalism both in the free case (on the space of infinite 
jets) and with constraints (on a PDE) are discussed. Analogs of tangent and 
cotangent bundles to a differential equation are introduced and the variational 
Schouten bracket is defined. General theoretical constructions are illustrated 
by a series of examples. 



Introduction 

The main task of this paper is to overview a series of our results achieved recently 
in understanding integrability properties of partial differential equations (PDEs) 
arising in mathematical physics and geometry |511 EOl EB E21 SH 113 ESI EH ES] ■ 
These results are essentially based on the geometrical approach to PDEs devel- 
oped since the 1970s by A. Vinogradov and his school (see [701 US [71 \IM W\ 
and references therein). The approach treats a PDE as an (infinite-dimensional) 
submanifold in the space J°°{t!') of infinite jets for a bundle it: E ^ M whose 
sections play the role of unknown functions (fields). This attitude allowed to apply 
to PDEs powerful techniques of differential geometry and homological algebra. The 
latter, in particular, made it possible to give an invariant and efficient formulation 
of higher-order Lagrangian formalism with constraints and calculus of variations 
(see [1301 [Hi [HI mZ]), see also [m d dUl [HI HH [751 [HS] and, as it became 
clear later, was a bridge to BRST cohomology in gauge theories, anti-field formalism 
and related topics, [IS]; see also [51 1151115] . 

Geometrical treatment of differential equations has a long history and originates 
in the works by Sophus Lie ^79,t80,.8L, as well as in research by A.V. Backlund [5], 
G. Monge [ST], G. Darboux [21], L. Bianchi [TO] and, later, by Elie Cartan [TO] . 
Note incidentally that Cartan's theory of involutivity for external differential sys- 
tems was an inspiration for another cohomological theory associated to PDEs and 
developed in papers by D. Spencer and his school, |115l 143) . Spencer's work (the 
so-called formal theory) closely relates to earlier and unfairly forgotten results by 
M. Janet [SB] and Ch. Riquier [TOT]; see also ^104] as weU as [TOBl [TTOl [HI [75] . 

A milestone in the geometry of differential equations was introduction by Charles 
Ehresmann the notion of jet bundles [331 135 that became the most adequate lan- 
guage for Lie's theory, but a real revival of the latter came with the works by 
L.V. Ovsyannikov (see his book |103j on group analysis of PDEs; see also [51] 11161 

mni]). 

A new impulse for the reappraisal of the Sophus Lie heritage was given by the 
discovery of integrability phenomena in nonlinear systems [3S1 US] UHl 11171 [351 
[Tllini [m ES] in the faff of the 1960fl3 and Hamihonian interpretation of this 



integrability [1381 [38] . In particular, it became clear that integrable equations 
possess infinite series of higher, or generalised, symmetries (see [1011 [70]). and 



"'^ Though discussions on "what is integrability" continued |137| and are still held now (see, e.g. 
quite recent papers I36| or I100| 'l. 
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classification of evolution equations with respect to this property allowed to discover 
new, at that time, integrable equations, [Ml 11181 155] . Later the notion of a higher 
symmetry was generalised further to that of a nonlocal one |133) and the search 
for a geometrical background of nonlocality led to the concept of a differential 
covering [TTj . The latter proved to play an important role in the geometry of PDEs 
and we discuss it in our review. 

It also became clear that the majority of integrable evolutionary systems possess 
a bi-Hamiltonian structure [531 [51], i-G-, can be represented as Hamiltonian flows 
on the space of infinite jets in at least two different ways and the corresponding 
Hamiltonian structures are compatible. The bi-hamiltonian property, by Magri's 
scheme [82| . leads to the existence of infinite series of commuting symmetries and 
conservation laws. In addition, it gives rise to a recursion operator for higher sym- 
metries that is an efficient tool for practical construction of symmetry hierarchies. 
Nevertheless, recursion operators exist for equations possessing no Hamiltonian 
structure at all (e.g., for the Burgers equation). A self-contained cohomological ap- 
proach to recursion operators based on Nijenhuis brackets and related to the theory 
of deformations for PDE structures is exposed in [51] . 

The literature on the Hamiltonian theory of PDEs is vast and we confine our- 
selves here to the key references [35l [25l [28l [30] , but one feature is common to all 
research: theories and techniques are applicable to evolution equations only. Then 
a natural question arises: what to do if the equation at hand is not represented in 
the evolutionary form? We believe that (at least, a partial) answer to this question 
can be found in this paper. 

Of course, one of possible solutions is to transform the equation to the evolu- 
tionary form. But: 

• Not all equations can be rendered to this forrr0. 

• How to check independence of Hamiltonian (and other) structures on a 
particular representation of our equation? In other words, if we found a 
Hamiltonian operator in one representation what guarantees that it survives 
when the representation is changed? 

• Even if the answer to the previous question is positive, how to transform 
the results when passing to the initial form of the equation? 

In what follows, we treat any concrete equation "as is" and try to uncover those 
objects and constructions that are naturally associated to this equation. In par- 
ticular, we do not assume existence of any additional structures that enrich the 
equation. Such structures are by all means extremely interesting and lead to very 
nontrivial classes of equations (e.g., equations of hydrodynamical type |3J, Monge- 
Ampere equations |77) or equations associated to Lie groups 29[), but here we look 
for internal properties of an arbitrary PDE. 

As it was said in the very beginning, an equation (or, to be more precise, its infi- 
nite prolongation, i.e., equation itself together with all its differential consequences) 
is a submanifold in a jet space J°°{n). To escape technical difhculties, we consider 
the simplest case, when tt: E ^ M is a. locally trivial vector bundle, though all the 
results remain valid in a more complicated situation (e.g., for jets of submanifolds). 
The reader who is interested in local results only may keep in mind the trivial 
bundle x R" M" instead of tt. 

Understood in such a way, any equation £ C J°°(7r) is naturally endowed with 
a (dimM)-dimensional integrable distribution C (the Cartan distribution) which 
consists, informally speaking, of planes tangent to formal solutions of £. This is 
the main and essentially the only geometric structure that we use. 



For example, gauge invariant equations, such the Yang-Mills, Maxwell, Einstein equation, 
etc., can not be presented in the evolutionary form. 
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In the research of the PDE differential geometry we use somewhat informal but 
quite productive guidelines which were originally introduced in [129) (see also [1311 
11281 ') and may be formulated as 

The structural principle: Any construction and concept must take into 

account the Cartan distribution on £. 
The correspondence principle: "Physical dimension" of f is n = dimC = 

dim M and differential geometry of £ reduces to the finite-dimensional one 

when passing to the "classical limit n ^ 0" . 
The invariance principle: All constructions must be independent of the 

embedding £ J°°{tt) and defined by the equation £ itself. 

Below we accompany our exposition by toy dictionaries that illustrate the corre- 
spondence between two languages, those of the geometry of PDEs and classical 
differential geometry. 

The paper consists of three sections. In Section [T] we describe the geometry of 
the "empty equation", i.e., of the jet space J°°(7r). In particular, we define the 
tangent and cotangent bundles to J°°{tt), introduce variational difi^erential forms 
and multivectors and define the variational Schouten bracket. We discuss geometry 
of Hamiltonian flows on the space of infinite jets (i.e., Hamiltonian evolutionary 
equations) and Lagrangian formalism without constraints. Section [5] deals with the 
same matters, but in the context of a differential equation £ C J°°{tt). Although 
the exposition in this part is quite general, the result on the Hamiltonian theory 
(the definition of the cotangent bundle, in particular) are valid for the so-called 
2-line equations only (we call such equations normal in Section [2]). This notion is 
related to the cohomological length of the compatibility complex for the lineariza- 
tion operator of £ and manifests itself, for example, in the number of nontrivial lines 
in Vinogradov's C-spectral sequence, [HO] (see also [HTl [HI [H] and [H [TM] ). 
From this point of view, jet spaces are 1-line equations and this is the reason why 
they have to be treated separately. Finally, in Section |3] we briefly overview the 
theory of differential coverings for PDEs and some of its applications: nonlocal 
symmetries and shadows, Backlund transformations, etc. 

Our exposition of general facts related to the geometry of jet spaces and infln- 
itely prolonged equations, including the nonlocal theory, is essentially based on the 
books [701 HSl El] ■ Lagrangian formalism, both in the free case (on jets) and with 
constraints (on equations), is exposed using the material of |132] (see also [Ml [72]). 
The geometrical approach to Hamiltonian formalism (including the theory of the 
Schouten bracket) is based on [5^IM] . 

A practical implementation of the general theory, in the majority of cases, needs 
the use of an appropriate computer algebra software. To avoid technical details that 
obscure the essentials, we chose for a "tutorial example" the well known Korteweg- 
de Vries equation for which all computations are transparent and can be done by 
hand. We did our best to illustrate the theory by a reasonable number of less trivial 
examples and really do hope that the result will be interesting to the readers. 

1. Jet spaces 

Jet spaces constitute a natural geometric environment for differential equations 
and for equations of mathematical physics, in particular. But these spaces are 
themselves an interesting geometric object that contains information on Lagrangian 
and Hamiltonian formalisms without constraints. Thus, we begin our exposition 
with a description of these spaces and structures related to them. 
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1.1. Definition of jet spaces. Let ir: E ^ M he a locally trivial smooth vector 
bundlcQ over a smooth manifold M, dimAf = n, dimi? = m + n. In what follows, 
M will be the manifold of independent variables while sections of tt will play the 
role of unknown functions (fields) . The set of all sections s: M ^ E will be denoted 
by r(7r) and it forms a module over the algebra C°°(Af). Two sections s, s' G r(7r) 
are said to be k-equivalent at a point x G M if their graphs are tangent to each 
other with order k at the point s{x) = s'{x) G E. The equivalence class of s with 
respect to this relation is denoted by [s]^ and is called the k-jet of s at x. The set 

= I xeM, serin)} 

is endowed with a natural structure of a smooth manifold; the latter is called the 
manifold of k- jets of sections of tt. Moreover, the maps 

TTfclj'^W^Af, (1) 

and 

TTfe,: J'=(7r)^ J'(7r), [sf, ^ [s]l k>l, (2) 

are smooth fibre bundles, tt^ being vector bundles. For any section s G r(7r) the 
map 

jfc(s): Af J'^'W, x^[st, (3) 

is a smooth section of tt^ that is called the k-jet of s. 

Here we are mostly interested in the case k = oo, i.e., in the space J°°(7r). It 
can be understood as the inverse limit of the chain 

. . . ^ IJ1±1^ jfc(^) ^ . . . ^ ji(^) :!i£, jO(^) ^ ^ M. (4) 

Due to projections Q there exist monomorphisms of function algebras 

C°^{M) c Fo{n) c • • • C Fk{n) c Fk+i{n) c . . . , (5) 

where J^fc(7r) — [J^ [tt)) , and we define the algebra of smooth functions on J°°(7r) 
as the filtered algebra J^(7r) = UfcJ^fc(7r). Elements of J^(7r) are identified with 
nonlinear scalar differential operators acting on sections of tt by the following rule: 

A/(s)=j^(s)*(/), sGr(7r),/G-FW. (6) 

More general, let tt': £" — > M be another vector bundle and 7r*(7r') be its pull- 
back to J°°(7r). Introduce the notation F{tt, tt') = r(7r*(7r')). Then any section / G 
J-'(7r, tt') is identified, by a formula similar to (jB]), with a nonlinear differential 
operator that acts from r(7r) to r(7r'). 

1.2. Vector fields and differential forms. A vector field on J°°{Tr) is a deriva- 
tion of the function algebra J-{n), i.e., an R-linear map X: TItt) — ^ F{tt) such 
that 

^(/5) = /^(5) + .9^(/) 

for all /, 5 G F{tt). The set of all vector fields is denoted by X{tt) and it is a Lie 
algebra with respect to the commutator (the Lie bracket). 

The definition of a differential form of degree r on J°° (tt) is similar to that of 
smooth functions. Using projections Q we consider the embeddings A''(j'^(7r)) C 
^r^jfc+i^^-j-j and set A'"(7r) = UfcA''( j'^(7r)). We shall also consider the Grassmann 
algebra of all forms A*(7r) — ©r>oA'"(7r) with respect to the wedge product. 



'For a definition of jets in a more general setting see, e.g., 1701 11321 IllOl [TTT) . 
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Coordinates. Let U C M he a coordinate neighbourhood such that the bundle tt 
becomes trivial over U. Choose local coordinates x^, . . . ,x" in U and u^, . . . , 
along the fibres of tt over U. Then the adapted coordinates in ■k~^(U) C J°°{'k) 
naturally arise. These coordinates are denoted by Mj, / being a multi-index, and 
are defined by 

j«.(s)*K) = 

where s = (s^, . . . , s™) is a local section of tt over U. In other words, the coordinate 
functions u'j correspond to partial derivatives of local sections. 
In these coordinates, smooth function on J°°{M) are of the form 

where the number of arguments is finite. Vector fields are represented as infinite 
sums 



i 1,3 I 

while differential forms of degree r are finite sums 

cj = V , dx'^ A • • • A Ax'" A Au]-'^' A • • • A duj" . 

1.3. Main structure: the Cartan distribution. Let 9 G J°°{tt). Then the 
graphs of all sections joo(s), s £ r(7r), passing through the point have a common 
n-dimensional tangent plane Cg (the Cartan plane). The correspondence C: ^ Co 
is an integrabl^l n-dimensional distribution on J°° (tt) that is called the Cartan dis- 
tribution. This distribution is the basic geometric structure on the manifold J°° (tt) . 
In particular, the following result is valid: 

Proposition 1.1. A submanifold in J°°{M) is a maximal integral manifold of C if 
and only if it is the graph of joo(s), where s is a local section of tt. 

Moreover, since the planes Ce project to M non-degenerately, any vector field X 
on M can be uniquely lifted up to a field CX on J°° (tt) . In such a way, one obtains 
a connection in the bundle tTqc called the Cartan connection. This connection is 
flat, i.e., 

C[X,Y] = [CX,CY] (7) 

for all vector fields X, Y on M. Due to (O, the space CX{tt) of all vector fields 
lying in the Cartan distribution is a Lie subalgebra in X{tt). Vector fields belonging 
to X{t:) will be called Cartan fields. 

Any vector field Z G X{Tr) can be uniquely decomposed to its vertical and 
horizontal components, 

Z = Z'' + Z^, (8) 

where Z'" is the projection of X to the fibre of the bundle Woo along Cartan planes, 
while Z^ lies in the Cartan distribution. Thus, one has 

X{tt) = X''{Tr)(BCX{Tr), (9) 

where ^'"{tt) is the Lie algebra of vertical vector fields. 

Dually to ©, the module of differential forms A^(7r) splits into the direct sum 

AiW=Ai(7r)©AiW, (10) 



^Integrability is understood formally here and means that if X and Y are two vector fields lying 
in C then their bracket [X, y] lies in C as well. Since J°°{tt) is infinite-dimensional, this does not 
mean that the Frobenius theorem holds for C: for any point 9 £ J°°{tt) there exist infinitely many 
maximal integral manifolds that contain 0. On the other hand, if C is restricted to an equation 
(see below Section [2]|, there may exist no maximal integral manifold at all. 
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where Aj(7r) consists of l-forms that annihilate the Cartan distribution (they will 
be called Cartan forms, or higher contact forms), while elements of Ajj(7r) are 
horizontal forms. 

Coordinates. Choose an adapted coordinate system (a;%Wj) in J°°(7r). Then one 
has 

The fields Di are called total derivatives and they span the Cartan distribution. 
For a basis in the module Aj(7r) one can choose the forms 

uj^j ^ du\ ~^ui^Ax\ (12) 

i 

while horizontal forms are 

ui — ''^^a.iAx^ , Oi e J^(7r). (13) 

i 

Remark 1 . It should be noted that all results and constructions below are valid not 
for the entire jet space only, but for an arbitrary open domain in J°°(7r). Everywhere 
below, when speaking about J°°(7r), we actually mean an open domain. 

1.4. Evolutionary vector fields and linearizations. We shall now describe 
infinitesimal symmetries of the Cartan distribution on J°° (tt) . A vector field X 6 
X{'k) is a symmetry if [X, Z] € CX{tt) as soon as Z e CX{tt). The space CX{tt) is 
an ideal in the Lie algebra Xc{'k) of symmetries. Due to integrability of the Cartan 
distribution, any Z G CX(t:) is a symmetry, and we call such symmetries trivial. 
Thus, we introduce the Lie algebra of nontrivial symmetries as 

syniTT = Xc{-K)/eX{-K). 

By ^ , sym tt is identified with the vertical part of Xc (tt) . 

Take a vector field X G symTr and restrict it to the subalgebra J-q{t^) C ^(tt). 
Then this restriction can be identified with an element ipx G -^(tt, tt). For shortness, 
we shall use the notation J^(7r,7r) ~ x{'k). 

Theorem 1.2. Tiie correspondence X i-^ ipx defines a bijection between symTr 
and x{'k). 

The element (px is called the generating section of a symmetry X, while the 
symmetry corresponding to a section S x{'k) is called an evolutionary vector field 
and is denoted by E^. 

Theorem 11.21 allows to introduce an 7^(7r)-module structure into the Lie alge- 
bra sym TT by setting 

f ■ E^= Ef^. 

This multiplication differs from the usual multiplication of a vector field by functions 
and does not survive when passing from the space of jets to equations (see Section[5] 
below). 

On the other hand, the same theorem defines a Lie algebra structure in x{t:): 
the Jacobi bracket {ip, ip} is uniquely given by the equality 

E{^^^} = [Eip,E^]. (14) 

The Jacobi bracket can also be computed using the formula 



{ip,i:}^E^{-4,)-E^{^). 



(15) 
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Remark 2 . We pointed out in Footnote|3]that integrability of distributions on J°° (tt) 
differs from tire one on finite-dimensional manifolds. The same holds for integra- 
bility (i.e., existence of the corresponding one-parameter group of transformations) 
of vector fields. Generally speaking, a vector field on J°°{tt) is not integrable in 
this sense, but there exists an important class of vector fields that are integrable. 
Namely, if X is a field on J''{tt), k < oo, that preserves the Cartan distribution 
then it can be lifted in a natural way to j'^+^(7r), see [13]. The entire collection 
of such fields determines a vector field on J°°(7r), and any such a field possesses 
a one-parameter group of transformations. The infinitesimal version of the Lie- 
Backlund theorem states that all such fields are the lifts of arbitrary fields on J'^{tt) 
(when dimTT > 1) or of a contact vector field on J^(7r) (when dimTr = 1). A 
complete description of integrable vector fields on J°° (tt) can be found in [121 120] ■ 

Coordinates. Let, in an adapted coordinate system, a section if g ^{tt) be of the 
form (p = ((/3^, . . . , 93™). Then the corresponding evolutionary vector field is 

E^^J2^^(V')^^ (16) 

where Dj ~ Di-^ o • • • o Di^ is the composition of total derivatives corresponding to 
the multi- index I = ii . . .ii. 

li ip = {ip^ , . . . , Ip"'-) is another element of x{tt) then the components of the 
Jacobi bracket are 

i;y = Y: - Diir)^) , J - 1, . . . , m. (17) 

Fix a section ip G >c{tt) and consider the map 

: >c{tt) x{tt), i^{p>) = E^iiP). (18) 

This map is called the linearization of the element ip (recall that "0 may be identified 
with a nonlinear differential operator acting from tt to tt). 

More generally, let tt' : £" M be a vector bundle. Then the action of an 
evolutionary vector field can be extended to 

E^: J-(7r, tt') ^ ^"(71, tt') 

in a well defined way. Consider a section ip e J^(7r, tt'), i.e., a nonlinear differential 
operator from r(7r) to r(7r'). Its linearization is the map 

l^: x{tt) ^ T{TTy) (19) 

is defined similar to (IT51) . 

Coordinates. Let, in adapted coordinates, (p = {ip^ , . . . , ip'^) and ip = {ip^ , . . . , ip"^ ). 
Then the jth component of iif,{(p) is 

I, a ^ 

i.e., the linearization is a matrix operator of the form 

j — l,...,m' 



(20) 

a— 1, . . .,m. 

Using hnearizations, formula (jl5p can be rewritten as 

{^M = - W)- (21) 
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1.5. C-differential operators. From ([20]) we see that linearizations are differential 
operators in total derivatives. We call such operators C-differential operators. More 
precisely, let ^ and ^' be two vector bundles over J°°{'k) and P, P' be the J-"(7r)- 
modules of their sections. An R-linear map A : P — ^ P' is a C-differential operator 
of order k if for any point G J°°{tt) and a section p ^ P the value of A(p) at 6 is 
completely determined by the values of Di{p), \I\ < k, at this point. The space of 
all such operators is denoted by Cfe(P, P') and we also set C{P, P') — LlkCk{P, P')- 
A closely related notion to that of a C-differential operator is horizontal jets. 
Let P be as above. We say that two sections p, p' P are horizontally k-equivalent 
(the case fc = cx) is included) at a point 9 € J°°{tt) if Di{p) — Dj{p') at 9 for all / 
such that |/| < fc. Denote the equivalence class by {p}g. The set 

JHiP) = {{p}o I 9eJ°^in),peP} 

forms a smooth manifold which is fibred over J°°(7r), 

4^(P)^ J°°(7r) {pj'e^e. 

The section j^ip) : J°°(7r) — > J^(P), 9 {p}g, is called the horizontal jet of p € P. 

Proposition 1.3. For any C-differential operator A S Ck{P,P') there exists a 
unique morphism $a of vector bundles £,k and suci that A{p) = <&A(jfc(p)) for 
any p E P. 

Two natural identifications will be useful below. 

Proposition 1.4. For any vector bundle it one has: 

(1) The module Aj(7r) is isomorphic to C(>^(7r), J^(7r)). 

(2) The module J^{>i{TT)) is isomorphic to X^{tt). The vector fields corre- 
sponding to sections of the form jj^ {p) are evolutionary fields. 

Coordinates. Choose an adapted coordinate system in the manifold J°°(7r) and let r, 
r' be dimensions of the bundles ^, respectively. Then any operator A G C(P, P') 
is of the form 

0=1,. ..,r 

ai^ e Tin). (22) 

a— l,...,r'. 

If w^, . . . jW*" are fibre- wise coordinates in the bundle ^ then the adapted coor- 
dinates in Jj^{(,), K being a multi-index, I — 1, . . . ,r, are determined by the 
equalities 

jLisnv'K)^DKis'), (23) 
where s = (s^, . . . , s'') is a local section of the bundle ^. 

Remark 3. The space of horizontal jets J^{P) is also endowed with an integrable 
distribution similar to the Cartan one: if 6* G J^{P) then the corresponding 
plane Ce is tangent to the graphs of horizontal jets passing through this point. The 
differential of the map ^oo : Jh°{P) ~^ isomorphically projects Ce to C^^(^gy 

Moreover, if P is of the form P — r(7r^(^)), where ^ is a vector bundle over M 
and 7r^(^) is its pull-back, then one has a diffeomorphism 

J^riP)^J°^{7TXM0, 

where tt Xm C is the Whitney product, and this isomorphism takes the Cartan 
distribution on J}^ (P) to the one on J°° (tt Xm 0- 

Remark 4. In the case when the modules P — r(7r^(^)) and P' = r(7r^(^')) are 
of the form considered in Remark [3] C-differential operators from P to P' may 
be understood as non-linear differential operators that take sections of tt to linear 
differential operators acting from ^ to 



A = 



a/3 



Dr 



GEOMETRY OF JET SPACES AND INTEGRABLE SYSTEMS 



9 



1.6. The variational complex and Lagrangian formalism. Consider now de- 
composition (|10p and note that it implies a more general splitting 

A^'W^ AgW^A^W, (24) 

p+q=k 

where 

A^(7r) = AJ W A • . • A AJW, A^Jtt) = Ai(7r) A ■ • • A Ai{n) . 
^ ^ ^ ^ ^ ^ 

p times q times 

Introduce the notation A^(7r) ig) A^^(7r) = EI''^{-k). By Proposition [Til (1), this 
space is identified with the module Cp'^(>f(7r), A^(7r)) of p-linear skew-symmetric 
C-difFerential operators acting from >f(7r) to A^(7r). 

The de Rham differential d: K^{ti) — > A'^+^(7r), by (p4l) . splits into two parts d = 
dc + d?i , where 

de = d?'-^: £;o^'«W ^ i?o^+^^^(7r), d, = d^, ^ : £;oP'«W ^ (25) 

are the vertical (or Cartan) and horizontal differentials, respectively These differ- 
entials anti-commute, i.e., 

dc o d,, + d,, o dc = 0. (26) 

Coordinates. For a function on J°°(7r) (i.e., a 0-form) the action of the Cartan 
differential is given by 

where dcUj = are the Cartan forms presented in p2p . while the horizontal 
differential acts as follows 

d^/ = ^A(/)da:\ (28) 



To compute the action on arbitrary forms it suffices to use ([27]) and (j28p and the 
fact that dc and d;i differentiate the wedge product and anti-commute with the 
de Rham differential. 

Let i?f''(7r) be the cohomology of at the term Eq''^{it), i.e., 

El'\^) = kerdf;7imd^'«"\ (29) 

Due to the vertical differential dc induces the differentials i?f '"'(tt) E\'^^'''{'k) 
which will be denoted by 5^'''. The groups i?f''^(7r), together with the differen- 
tials (JP'*, play one of the most important roles in the geometry of jets providing 
the background for Lagrangian formalism without constraints. To describe them, 
we shall need new notions. 

Let ^ be a vector bundle over J°°{'k) and P — T{^). Introduce the adjoint 
module P = \lou^Jr^^■^{P^ K^{-k)) . Consider another module of sections Q and a C- 
differential operator A: P ^ Q. Then the adjoint operator A* : Q — >■ P is defined 
and it enjoys the Green formula 

{A{p),q)-{p,A*{q))=d^io (30) 

for all p e P, q & Q and some oj — io{p,q) G A]^^^{'k), where (• , ■} is the natural 
pairing between the module and its adjoint. It is useful to keep in mind that the 
correspondence {p,q) H> Lu(p,q) is a C-diffcrcntial operator with respect to both 
arguments. 

An operator A is called self-adjoint if A* = A and skew-adjoint if A* = —A. 
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Coordinates. If A = ajDi is a scalar operator then 

A* =^(-1)1^1/?/ o a/. (31) 

For a matrix C-difFerential operator A = ||A,j || one has 

A* - II A*, II , (32) 

where A*j is given by ([5T|) . 

We can now describe the groups i?f'''(7r). 

Theorem 1.5 (One-hne Theorem). Let tt: E ^ M, dim M = n, be a locally trivial 
vector bundle. Then: 

(1) the groups E°''^{Tr), q — 0, . . . ,n — 1, are isomorphic to the de Rham coho- 
mology groups H'^{M) of the manifold M; 

(2) the group E^'"{tt) consists of the Lagrangians depending on the fields that 
are sections of tt; 

(3) the groups E^'^ (n) , p > 0, are identified with the modules C^*i{}<{'k) , >c{t:)) 
of {p — l)-linear skew-symmetric C -differential operators that are skew- 
adjoint in each argument (in particular, eI'"'{it) — ^({tt)); 

(4) all other terms are trivial. 

Remark 5. The group _E^'"(7r) is also called the nth horizontal cohomology group 
of TT and denoted by HJ^{tt). 

Remark 6. The construction above is a particular case of Vinogradov's C-spectral 
sequence, [l26l [T27l [T30l [T32l [72] . 

In what follows, we shall assume the manifold M to be cohomologically trivial, 
i.e., its de Rham cohomology is isomorphic to E. 

Define the operator 6 : A5J(7r) — k{TT) as the composition of the projection A5J(7r) — 
and the differential 5?'": ^i"^" ^ 

Remark 7. In what follows, we shall also use the notation 6 for the differential Si'" 
itself. 

Coordinates. If w G A^(7r) is of the form L dx^ A • • • A dx" then 
where 

are variational derivatives. Thus, S is the Euler operator and it takes a Lagrangian 
density w to the corresponding Euler operator. 

Proposition 1.6. Let tt be a locally trivial vector bundle over a cohomologically 
trivial manifold M . Then the complex 

Tin) ^ AUt:) ^ . . . ^ Al-\n) A^(7r) 

^ k^Cf*{>i,k) ^C^^*{x,k) ^ ■■■ (35) 

is exact, i.e., the kernel of each differential coincides with the image of the preceding 
one. 

In the coordinate-free way, the Euler-Lagrange operator can be computed by 
S{u) =£:(!), ioeAU^), (36) 
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while the difFerentials S^'"' enjoy the equahty 

1=1 

(37) 

where A G C^*i{K{Tr), >c{Tr)) and ipi, . . . ,ipp G >c. Here and below we use the 
notation ^A,i^i,...,i^fc = E^{A){ipi, . . . , ipk) In particular, if ip G >ir(7r) then 

=^^-^;. (38) 

Remark 8. Complex (1351) is exact starting from the term A5J(7r) independently of 
cohomological properties of the manifold M. This complex is called the global 
variational complex of the bundle tt, see [134) . 

As a consequence of Proposition II .61 we obtain the following result: 

Theorem 1.7. For a vector bundle tt one has: 

(1) The action functional 

Jm 

is stationary on a section s if and only if joo(s)*((5(a;)) = (i.e., as we shall 
see below, s is a solution of the Euler-Lagrange equation corresponding 
to uj). 

(2) Variationally trivial Lagrangians are total divergences, which amounts to 
the equality keri5 = imd^'"~ . 

(3) All null total divergences are total curls, i.e., d°'"^^a; = if and only 
ifcu ^ d°'"~26l for some 9 e Al'^^Tr). 

(4) If ^ G k{TT) then the nonlinear differential operator : >c{tt) k{TT) is 
of the form 5u) (i.e., is an Euler-Lagrange operator) if and only if = £^ 
(the Helmholtz condition). 

1.7. A parallel with finite-dimensional differential geometry. I. We want 
to indicate here a very useful and productive analogy between the geometry of 
jet spaces (and, more generally, of differential equations) and classical differential 
geometry of finite-dimensional smooth manifolds. This parallel was exposed by 
A.M. Vinogradov first within his philosophy of Secondary Calculus (cf. [132j and 
references therein) and we elaborate it further. 

Two points of view on jet spaces as geometrical objects may exist. The first one 
is formal, traditional and straightforward. It was described on the previous pages 
and treats J°°(tt) as a particular case of general infinite-dimensional manifolds. 
Such an approach, being by all means necessary for a rigorous exposition of the 
theory, actually ignores essential intrinsic structures of jet spaces actually. 

Another viewpoint is completely informal but incorporates these structures ab 
ovo and allows to reveal new and non-trivial relations and results just "translat- 
ing" from the language of the classical differential geometry. To facilitate such 
translations, we shall compile a sort of a dictionary. 

So, we consider the space J°°(7r) endowed with the Cartan distribution C and 
take for the points of the new "manifold" maximal integral submanifolds of C As 
it was indicated above, they are graphs of joo(s), s € r(7r), and thus new "points" 
are sections of tt (i.e., fields). 

Let UJ S A5J(7r) be a horizontal n-form on J°°(7r) (or a Lagrangian density). Then 
to any "point" s G r(7r) we can put into correspondence the number 

W(s) = / joois)*(uj). 

Jm 
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Thus, Lagrangians are understood as functions. Due to the identity 

joo{s)*{dhUj) = d(joo(s)*(w)) 

and the Stokes formula 

[ d0= [ e, e e A"-i(Af), 

Jm JdM 

"functions" of the form d/iOj vanish at all points. So, no-trivial functions are el- 
ements of the cohomology group £^J'"(7r) = iJ^(7r). Thus, the beginning of the 
dictionary is 

Manifold M Jet space J°°(7r) 

points < — > sections of tt (fields) 
functions < — > Lagrangians oo — L dx^ A • • • A dx" 



value at a point, f{x) < — > integral uj{s) ~ / ioo(s)*w 

Jm 

(the cohomology class of uj) 

The next step is to define vector fields. They should be infinitesimal transfor- 
mations of J°°{tt) that preserve the Cartan distribution (or, equivalently, move 
"points" to "points"). These are exactly the fields lying in Xc{Tr). But vector 
fields X G CX{'!t) (i.e., lying in the Cartan distribution) are tangent to maximal 
integral manifolds of the latter and thus are trivial in the space of fields. Con- 
sequently, non-trivial vector fields are identified with elements of symTr, i.e., with 
evolutionary vector fields on J°°{tt). On the other hand, as it was indicated above, 
they are integrable sectionf0 of J^(x). Hence, this bundle can be considered as 
the tangent bundle to J°°(7r). So, the dictionary can be continued as follows: 

Manifold M Jet space J°°(7r) 

vector fields, X{M) i — > evolutionary vector fields, m:{tt) 
the tangent bundle i — > the bundle of horizontal jets J^{k{tt)) 

Differential forms on a smooth manifold M may be understood as multi-linear 
functions on the space of vector fields (or fibre- wise multi-linear functions on T{M)): 
we insert a vector field into a p-form and obtain a (p — l)-form. In our context, 
such objects are exactly elements of C^_*^{h, k) = £'f '"(tt). We call them variational 
forms of degree p and have the following parallel: 

Manifold M Jet space J°°(7r) 

differential forms, AP(Af) i — > variational forms, Cp'l*(>f, iir) 
the de Rham complex < — > the variational complex 

Remark 9. It can also be shown that smooth maps J°°(7r) — > J°°{tt') preserving 
the Cartan distributions are completely determined by non-linear differential op- 
erators from r(7r) to r(7r') while the differentials of these maps Jf^{K) Jf^{x') 
correspond to linearizations. Unfortunately, a detailed exposition of this parallel is 
out of scope of our review. 

We shall continue to compile our dictionary in Subsection 11.91 



By integrable sections we mean those ones whose graphs are maximal integral manifolds of 
the Cartan distribution. 
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1.8. Hamiltonian formalism. The objects dual to Cp'i2(>!r(7r), i<(7r)) are the mod- 
ules of variational multivectors Dp{Tr) = C^*i{>c{tt) , >c{tt)) . In particular, Di{tt) = 
>c{it). We also set Dq{'k) = A'^{Tr)/ imd^'"^^ — £'"'"(7r). To describe Hamiltonian 
formalism on J°°(7r), we first introduce the variational Schouten bracket j59j 

in the following way (cf. [68l[55], see also [281111]). If i? = [uj] £ Do{Tr) is a coset of 
a horizontal form u G AJJ(7r) and A G Dp{Tr), p > 0, then we set 

lA,Bj = {-inB,Aj^AiSu:), 

while for any B £ Dq{Tr), q > 0, and i/; = Slu <E >c{Tr), o-i G A^{tt), 

[A^IW = lA,Bii;)j + {-ir-'lA{i;),Bl 

and these two equalities define the bracket completely. In particular, 

and 

Proposition 1.8. The variational Schouten bracket determines a super Lie algebra 
structure in the space D^n) — J2p>o ^p(^) ^'^ ^^'^ following sense: 

= (39) 
(_l)(p+i)(-+i) 11^^ B], C] + (-l)('?+i)(p+i) IIB, ClAj (40) 
+ (_l)(-+i)(g+i)||C,Al,Bl =0 
for all A e Dp{Tr), B e Dg{Tr), C £ Dr{n). 

To compute the Schouten bracket explicitly, for any natural n consider the set 5*,^ 
of all (n — i)-un-shujfles consisting of all permutations a of the set {1, . . . , n} such 
that 

cr(l) < • • • < <7{i), a{i + 1) < ■ ■ ■ < a{n). 
We formally set 5,\ = for i < and i > n. We also use a short notation i^aikiM) 

for 1pa{ki),- ■ • ,'0<T(fc2)- 

Let now A e Dpi?:) and _B G Dq{T:). Then for any -01, . . . G >f(7r) we 

have 

|A, 51(01,..., Vp+g-l) - J2 (-ir^S,^.(i„_i,(A(0,(,,p+g-l))) 

_ ^ (-l)-i?(^^^^^^^^_^^(0,(,)),0,(,+ l,,+,_l)) 

_(_1)(P-1)(.-1) ^ (-l)"^A,^.<,.,_„(i3(0.(p.p+,-l))) 

+ ^ (-l)M(£^,,^^^^^^_^,(V-.(,)),0.(,+i,,+,_i)), (41) 

p + q-l 

where (—1)*^ stands for the parity of the permutation a and, as before, iA,ipt,...,ipk iv) — 
i;^(A)(0i,...,0fc). 

We say that a bivector A G D2{tt) = C^^*{>c{tt), >c{tt)) is a Hamiltonian structure 
on J°°(7r) if 

IAA1=0. (42) 

Remark 10. A more appropriate name is a Poisson structure but we follow here 
the tradition accepted in the theory of integrable systems. 
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Given a Hamiltonian structure, one can define a Poisson bracket (with respect 
to the Hamiltonian structure A) on the set of Lagrangians: 

{ij,Lu'}A^ {A{S{uj))J{uj')), u,Lu' e Do{n). (43) 

Two elements are in involution (with respect to the structure A) if 

{uj,ll>'}a = 0. 

Proposition 1.9. For any A S C^^* {>c{tt), >c{tt)) one has 

{uj,uj'}a^ -W,uj}a (44) 
If in addition A satisfies (|42|) then 

{uj, {uj', lu"}a}a + W, {u",u}a}a + W\ {lo, lo'}a}a = 0. (45) 
One also has 

A{u^u'}^ ^ {A^,A^,}, (46) 

where A^ = A(6(ijj)) G >f(7r) and the curlies in the right-hand side denote the 
Jacobi bracket. 

Chose a Hamiltonian structure A and consider the sequence of operators 

^ Doin) ^ Di{^) ^ > Dp{n) ^ Dp+iin) ^ • • • , (47) 

where 9^(5) = 

Proposition 1.10. Sequence (|T7)) is a complex, i.e., Oa o Qa = 0. 

We say that is a Hamiltonian vector field if e herd a- This is equivalent to 
E^{{lj,lu'}a) = {E^iuj),uj'}A + {lo, E^{uj')}a, (48) 

i.e., Eip preserves the Poisson bracket. Due to Proposition [TTTOl a particular case of 
Hamiltonian fields are fields of the form Ea(s(uj)}- In this case, ut G Do{tt) is called 
the Hamiltonian of the field under consideration. 

We say that oj e Do{Tr) is a first integral of a Hamiltonian field E^ if Eip{uj) = 0. 
A Hamiltonian field E^i is a symmetry for the field E^p if [E^pr, E^] — 0, or 

{E^-l^W)^Q. 

Proposition 1.11. If E^p is a Hamiltonian vector Geld with respect to a Hamilton- 
ian structure A then the operator A o S ^ Oa : -Do (ti") ^ -Di (tt) takes Erst integrals 
of E^ to its symmetries. 

A Hamiltonian structure B G D2 (tt) is said to be compatible with the structure A 
if B G ker Oa , or 

[A 51=0. (49) 

This is equivalent to the fact that all the bivectors 

AA + /iS, A,^gR, (50) 

are Hamiltonian structures on J°°{it). The family (jSOp is called a Poisson pencil. 
When two Hamiltonian structures are given, one also says that they form a bi- 
Hamiltonian structure. 

Coordinates. Let us indicate how to verify conditions (H^ and in coordinates 
(the explanation will be given below in Subsection II. 9[ see Remark [Tl)) . Take 
bivectors A, B £ D2{tt) — Cf^* {k{n), >c{tt)). Then A and B, in adapted coordinates 
in J°°{tt), are represented as matrix C-differential operators 



B 
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where i, j = 1, . . . ,m = dimTr. Let us put into correspondence to these operators 
the functions 

Wa = ' = E ' (51) 

where are odd variables. Then ^ is a Hamiltonian structure if and only if 



5Wa SWa 



while two Hamiltonian structures are compatible if and only if 



0, (52) 



'v- ( SWa 



6Wb 6Wb SWa\\ _ 



\ Su^ Sp^ Su^ 6p^ J 



0. (53) 



Theorem 1.12 (the Magri Scheme, see [HIEH])- Let {A,B) be a bi-Hamiltonian 
structure on J°°{tt) and assume that the complex ((47|l is acychc in the term Di{tt), 
i.e., every Hamiltonian vector Geld with respect to A possesses a Hamiltonian. 
Assume also that two densities wi, W2 G Do{Tr) are given, such that dA{uJi) — 
93(^2) ■ Then: 

(1) There exist elements W3, . . . , Wg, • • • G Do{tt) satisfying 

dAiuJs) = dsii^.+i), s = 2,3,... (54) 

(2) All elements uJi, . . . ,ujs, . . . are in involution with respect to both Hamil- 
tonian structures, i.e., 

{uJa,U!i3}A = {uJa,UJp}B = 

for aU a, I3>1. 

Example 1 (the KdV hierarchy). Consider J°°(7r) for the trivial one-dimensional 
bundle tt : R x K — M. Let x be the independent variable and u be the fibre coordi- 
nate (the unknown function). Then the adapted coordinates u — Uf), ui, . . . ,Uk, . ■ . 
in J°°{tt) arise, where Uk corresponds to d^u/dx^. The operators 

A r. d ^ d 



and 



dx f^^ duk 



B ^ Dl+4uD.^ + 2ui, 



as it can be easily checked using (j52p and (j53p . constitute a bi-Hamiltonian structure 
on J°°(7r). Then obviously for the horizontal forms 

wi — —u^dx, UJ2 — —udx 
2 2 

one has dAt^i — A{u) — ui and 83(^2 = B{^) = ui, i.e., 

OaUJi = dBUJ2- 

The first cohomology group of Oa is trivial (see 2J), and consequently we obtain 
an infinite series of first integrals and the corresponding symmetries. The second, 
after ui, symmetry is 6uui + u^: 

6uui + — OaHu^ — 2^1) ^2^) = 9b{—u^ dx). 
The corresponding flow on J°° (tt) is governed by the evolution equation 

Ut = QuUx + Uxxx': 
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thus, we obtain the Korteweg-de Vries equation and the corresponding hierarchy 
of commuting flows (the higher KdV equations). The entire family of commuting 
flows can be obtained by applying the Lenard recursion operator (see [40p 

R = B o A-^ = Dl+4u + 2uiD-^ (55) 

to the right-hand side of the first flow ut — sufficiently many times. 

Example 2 (the Boussinesq hierarchy). Consider the adapted coordinates x, u, 
V,. . . ,Uk, Vk, - ■ ■ in the space J°°(7r), where tt: x R — R is the trivial two- 
dimensional bundle over R. Then the operators 

, f DA „ _ faDl + uD, + iui ^vD, 



where cr is a real constant, form a bi-Hamiltonian structure. For the differential 
forms uJi = 2(u -|- w) da; and uj2 = uvdx one obviously has Oa^^i = 831^2- The 
arising hierarchy of commuting flows corresponds to the evolution equation 

Vt ^ Ux+ VVx 

which is the two- component Boussinesq system which can be obtained from the 
Kaup equation, see [SS]. Note that there exists another Hamiltonian operator 

where 







3 

1- l^^^viDl + {av2 


+ uv)Dx + ^(o'Ws -f uvi + uiv), 


^uv 






1 

2 ' 


^vu 




h {u+]v^)D^ + 
4 


-{Ui +OTi), 






1 






which is compatible both with A and B. In this sense, system (j56l) is tri- Hamiltonian. 

Example 3 (the KdV hierarchy, II). Let tt: R^ x R^ -j- R^ be the trivial three- 
dimensional bundle with the coordinates t in the base and u, v, w in the fibre. 
Introduce the adapted coordinates Ufe, Vk, Wk, where fc = 0, 1, 2, . . . , in J°°(7r) and 
consider the operators 

-2u -Dt - 2v 

A = I I -6m I , B = I 2u Dt -12u2 _ 2w 

'Dt + 2v 12^2 + 2w 8uDt + 4ui _ 

which form a bi-Hamiltonian structure on J°° (tt) . It can be easily seen that Oai-Oi = 
83^2, where 

1 3 1 

oji — [uw — -u^ -I- 2u^) di, UJ2 — — ( 2^^^ + 

Thus, we obtain a hierarchy of commuting flows whose first term is 

Ux = V, Vx ^ W, Wx = Ut - 6uv, 

which is obviously the KdV equation rewritten in a different way (cf. the pa- 
per inn])- 
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Note that the Lenard recursion operator (|55p in the new representation of the 
KdV hierarchy acquires the form 

/ -2u ~Dt - 2v 

R = i 2u Dt -12m2 _ 2w 

\-Dt + 2v Uu'^ - 2w 8uDt + 4ui 

Remark 11. In a recently pubhshed paper [6], the authors formulate a much weaker 
than triviality of the first Poisson cohomology group criterion for feasibility of the 
Magri scheme. The criterion is given in the framework of Dirac structures |28[ 123) 
(though it also admits a self-contained formulation) that unfortunately lie beyond 
the scope of our review. Note that Dirac structures that merge the notions of sym- 
plectic and Hamiltonian operators constitute an interesting object for geometrical 
research in PDEs. 

Remark 12. Normal forms for Hamiltonian operators of order < 5 and a "varia- 
tional" analog of the Darboux Lemma were presented in [31 13] . In |114| , one can 
find normal forms for operators of order < 11 and some classification results for 
operators of higher order. 

1.9. A parallel with finite-dimensional diflferential geometry. II. Let us 

now continue to construct the dictionary started in Section [1.71 

Of course, variational multivectors introduced in Section [TT51 are exact counter- 
parts of classical multivector fields in differential geometry. These fields are nat- 
urally understood as smooth functions on T*M skew-symmetric and multi- linear 
with respect to fibre variables. Exactly the same interpretation is valid for varia- 
tional vectors if one considers the bundle r* : J^{k) — )■ J°°(7rj^. Thus, we have 
the following translations: 

Manifold M Jet space J°°(7r) 

multivector fields < — > variational multivectors, Cp'l''i^(>!:, >r) 

Schoutcn bracket < — > variational Schouten bracket 

Poisson structure < — > Hamiltonian operator 

cotangent bundle, < — > variational cotangent bundle, 

T*M^M i — > J^{k) ^ J°°iTr) 

As it is well known, the tangent space T*M is endowed with the natural sym- 
plectic structure ft = dp A dq G A^(M) and, in addition, ft = dp, where the 
form p = pdq is defined invariantly as well. Similar constructions exist on J'f^{k). 
Let us show this. 

To this end, recall (see Remark [3]) that the manifold J^{k) is diffeomorphic 
to J°°(7r Xm '^), where tt = Hom(7r, /\" r*M). Hence, the module of variational 
1-forms on (tt) is isomorphic to 

k{TT Xm tt) = k{Tr) x joo(^) k{tt). (57) 

Then the 1-form (the analog of p dq) is uniquely defined by the condition 

jr(V')*(p.) = (^,0), (58) 

where "0 G k^n) is an arbitrary variational 1-form on J°°{t:). 

Now, by the same reasons and dually to ([57]) , the module of vector fields on (k) 

is 

M:{Tr Xm tt) = K{Tr) x,7oo(^) n{Tr). 



^Note that in such a way we independently arrived to Kupershmidt's notion of the cotangent 
bundle to a vector bundle, see 1751 . 
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Thus we see that the symplectic structure £7^ must be an element of the mod- 
ule C^^*{>c{tt Xm tt), >c{tt Xm tt)). For any element {ip,ip) G >c{tt Xm tt) we set 

n^iip,^) = i-^,ip); (59) 

this is a skew-adjoint C-differential operator of order 0. 

Remark 13. The form p^r can be defined in a different way. Namely, let X = {ip, tp) 
be a vector field on Jj^{k). Then its 1st component is a vector field vertical with 
respect to the projection J^{k) — > J°°{Tr) and may be understood as a function 
on J^{k). Then we set Ptt{X) = ip. This definition is equivalent to ((58)) . 

Of course, the operator fJ^r is invertible and the inverse one S^^ = Q,^^ is a 
bivector on the cotangent manifold J'f^{>c). There exist two points of view at 
this manifold (as well as at T*M). The first one treats it as a classical ("even") 
manifold. The second approach considers J^{>c) as a super-manifold with odd 
coordinates along the projection J^{k) — )• J°°(7r) and even ones in the base. If 
one takes the first approach the bivector St^ will define the Poisson bracket for 
functions on J^{>c). The second approach leads to functions multi-linear and skew- 
symmetric with respect to fibre variables. As it was stated above, these functions 
are identified with variational multivectors on J°°{tt). Then the super- bracket 
defined by S-^ coincides with the Schouten bracket. The bracket is given by the 
formula 

SASc.,)i.,) = {SASo^^lSc.,) = ( f'^^''^^/' 7; (60) 

1^ |[wi,W2j, the odd case, 

in both cases. 

Remark 14. The above said clarifies the meaning of formulas ([51 ]) - ([55]) . Namely, 
the correspondence A M- Wa given by ([5T|) describes how to construct the function 
on J^{k) when a bivector A is given (to be more precise, this function is the 
horizontal cohomology class of the form Wa da;^ A - • - Ada;" in HJ^{7t)). The argument 
of 6 in ([55)1 is the coordinate expression of the bracket (ISO)) in the odd case while (|Sn|) 
itself checks triviality of its horizontal cohomology class. 

2. Differential Equations 

With the concept of the jet bundle at our disposal we give a geometric definition 
of differential equations. 

2.1. Definition of diflferential equations. Suppose we have a system 

F,{x\u^)^0, s = l,...,/, (61) 

of partial differential equations in n independent variables and m dependent 
variables . Equations (pT|) determine a locus in the jet space J°° (tt) of a vector 
bundle tt: E ^ M, such that dimi? = m + n, dimM = n. 

The subset of (tt) defined in this way is not an adequate geometric construc- 
tion corresponding to the system at hand, because it does not take into account 
differential consequences of (|ST]) . So, we extend ([CT|) to a larger system 

Di{Fs) = for all multi- indices / and s = 1, . . . ,1, (62) 

and consider the locus £ C J°°{t!') defined by (|62p. 
Thus, we get a correspondence 

This correspondence behaves nice with respect to solutions of (|6ip : they are those 
sections of tt whose infinite jets lie in £. To put this another way, the solutions 
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of (pT|) are the maximal integral submanifolds of the Cartan distribution restricted 
to £. 

This shows that £ endowed with the Cartan distribution can be taken as the 
geometric object corresponding to system (j61l) . we call such a manifold £ an equa- 
tion. 

An equation is generally of infinite dimension. 

Without loss of generality we assume that (IHlT) does not contain equations of zero 
order, in geometric language this means that the projection tTocoI^ '■ £ ~^ J^i"^) is 
surjective. It is obvious that at every point 6 d £ the Cartan plane is tangent to 
the equation, Cg C Tg{£), so that the dimension of the Cartan distribution on an 
equation is equal to n, the same as on the jet space. 

System (1611) is a coordinate description of an equation £. Every equation has 
many different coordinate descriptions. 

Example 4. Take the KdV equation (cf. Example [T] from Section [1]) 

Ut - Guuoo - Uxxx = 0. (63) 

The bundle tt here is the projection tt: M'' — > R^, {x,t,u) H> {x,t). The jet space 
J°°(7r) has coordinates x, t, u, Ux, ut, . . . , uj, ... The equation £ C J°°(7r) is 
given by the infinite series of equations 

Ut = 6uux + U.^xx, 
UtI = Di{6uUx + Uxxx), 



The Cartan distribution on £ is two-dimensional and generated by 

s 

d d 
A = ^ + ^ Dli&uux + u^xx)——, 

s * 

where u^ = Ux...x {s times). The functions x, Us can be taken to be coordinates 
on £. 

The system 

Ux — V — 0, Vx — w ~ Q, Wx ~ Ut + Guv ~ (64) 

gives rise to the same equation £ C J°°(7r') with tt': M^, {x,t,u,v,w) 

{x,t). To prove that and (IMl) define the same equation consider the map 

a: J^(7r) ^ J°{tt'), a{x,t,u,Ux,ut,Uxx,Uxt,Utt) = ix,t,u,Ux,Uxx), 

b: J°{tt') J"(vr), b{x,t,u,v,w) = {x,t,u). 

Let a°°: ■/°°(7r) ^ J°°(7r') and 6°°: J°^{tt') J°°{n) be the lifts of these maps 
(cf. Remark [5]). Then it is easy to see that a°°|^o 6°°|g = b°°\£° 0'°°\s — i^l^- Lifts 
preserve the Cartan distributions, hence the maps a°°\g and b°°\^ are isomorphisms 
of equations determined by ([55]) and 



20 



JOSEPH KRASIL'SHCHIK AND ALEXANDER VERBOVETSKY 



Thus, two different coordinate expressions and ([M)) determine the same 
equation £ included into two different jet spaces: 

J°°(7r) 



£ 




J°°{tt'). 

As usual, the choice of functions (1611) should be restricted by some regularity 
assumptions. Namely, we require that there exists a subset S C {Di{Fs)} of 
functions on J°°{tt) such that 

(1) Fs e S for all s = 1, . . . , 

(2) the functions that belong to S define the equation £; 

(3) the differentials d/ are locally linearly independent on £ for all / G S. 

We always require these conditions to be satisfied. 

In this review, we shall assume that an equation at hand £ C J°°(7r) is globally 
defined by a relation F = 0, where F is a section of an appropriate /-dimensional 
vector bundle ^ over the jet space J°°{-k). This is always possible. Equations (|6T|) 
are local coordinate expressions for F — Q. Denote by P the J^(7r)-module of 
sections of ^, so that F G P. 

The above regularity conditions imply the following very useful fact: a function 
/ £ -/""(tt) vanishes on £, f\g = 0, if and only if / = A(_F) for some C-differential 
operator A: P — >■ J-"(7r). 

Let £ = {F = 0} be the equation defined by a section F ^ P. The section F is 
called normal if for any C-differential operator A: P — > -^('i') such that A(P) = 
we have A|g — 0. The equation £ C J°°(7r) is called normal if it can be defined by 
a normal section. 

Example 5. A simple example of an abnormal equation is the system 

Uy -V^ = 0, Uz -Wx = 0, Vz - Wy = 0. 

Gauge equations, including Maxwell, Yang-Mills, and Einstein equations, are not 
normal as well. Such equations are beyond the scope of our review. On the other 
hand, the majority of equations of mathematical physics, in particular all evolution 
equations, are normal. 

The Cartan vector fields on an equation £ form a Lie algebra CX{£). In the 
same manner as for jet spaces we define the Lie algebra of symmetries of £ and 
spaces of the Cartan and horizontal forms: 

syin£ = Xc{£)/CX{£), 

where Xc{£) = {X X{E) \ [X,CX{£)] C CX{£) }, 

AP{£) = {uj£ AP{£) I ixiij) =0 VX G CXi£) }, 

Al{£)^A\£)/Al{£), 

Al{£)^Al{£)A---AAl{£). 

We shall discuss them in more detail below. 

A morphism of equations f : £ ^ £' \s a smooth map that respects the Cartan 
distribution, i.e., for all points 9 ^ £ we have f*{Ce) C Cf(^g)^ where Cg is the Cartan 
plane at a point 6 ^ £. 

If a morphism f : £ ^ £' is a, fibre bundle and the map f*: Cg — ^ ^^(e) is an 
isomorphism of vector spaces for all points 9 G £ then / is called a covering. 

We shall discuss the theory of covering in Section [31 
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Example 6. For an equation £ we construct the quotient bundle 

r: T{£)^T{£)/C^£, 

where T{£) £ is the tangent bundle to f , C C T{£) is the Cartan distribution 
thought of as a subbundle of T{£). Given an inclusion £ C J°°(7r), tt: E ^ AI, 
as above, the fibre bundle r : T{£) — ^ £ can be identified with the vertical bundle 
with respect to the projection £ — > M. 

Every Cartan vector field X G CX{£) can be lifted to a vector field X e X{T{£)) 
as follows. It suffices to define an action of X on fibre- wise linear functions on T{£) 
that can be naturally identified with Cartan 1-forms uj G ^\{£)- We put 

where Lx denotes the Lie derivative. In coordinates, we have Di ~ Di, where Di 
in the right-hand side are the total derivatives on T{£). 

Let us describe an inclusion of T{£) to a jet space. Assume that £ C J°° (tt) is 
given by = 0. Let x- J°°{t^ x m tt) — > J°°(7r) be defined by the projection to the 
first factor. Then T{£) C J°°{Tr Xm tt) is defined by equations 

F = 0, iF{v)^0, 

where the tilde over F denotes the pullback by x (with the tilde for the Cartan 
vector fields defined above), v G >ir(7r) = T{n Xm 7'',7r) is the projection to the 
second factor r(7r x a/ vr) r(7r). 

So, T{£) is an equation and the vector fields of the form X generate the Cartan 
distribution on it. Thus, the bundle r : T{£) — >■ £ is a covering called the tangent 
covering to £. 

In coordinates, we have v — {v^, . . . , w™) if the coordinates on J°°{tt Xm tt) are 
x% Uj, Uj, with and v^'j corresponding to the first and second factors, respectively. 
Thus, a coordinate description of T{£) has the form 

a, I ^ 

Note that if f is a normal equation then T{£) is normal as well. 

The reader will find more examples and a detailed discussion of coverings in 
Section [3] below. 

2.2. Linearization. Let £ C J°°{tt) be an equation defined by a section F G P. 
Denote by k the restriction of the module >c('k) = J^(7r,7r) to £. The linearization 
of the equation £ is the restriction to £ of the linearization of F: 

Ip = lp\^ : m:^ P. 

We denote by bar the restriction of a C-difFerential operator to £ and preserve the 
notation of modules for their restrictions. 

Remark 15. The operator £p is well-defined globally only if the module P has the 
form P — 7^(7r,7r'). For an arbitrary module P the operator £p is defined only 
locally. But its restriction £p is well-defined globally on the whole £. 

Remark 16. A C-difFerential operator A on f is called normal if for any C-differential 
operator □ the condition □ o A = implies = 0. 

The linearization operator £p of an equation £ is normal if and only if the 
section F is normal. If an equation £ admits a normal representation, i.e., an 
embedding £ — {F = 0} C J°°(7r) such that the corresponding operator £p is 
normal then there always exists another representation £ = {F' = 0} C J°°{tt') for 
which £ acquires an evolutionary form. 
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Recall that in coordinates the linearization £p has the form (cf. ([^0])): 

j=i,...,i 



Q — l,...,m. 

If an equation S is defined by two different sections Fi G Pi an F2 G P2 in 
different, generally speaking, jet spaces 

J°°(^i) 

then the corresponding linearizations Ipi'- xi Pi and £p^: X2 ^ P2 are equivalent 
[311I53] in the sense that there exist C-differential operators a, /3, a', Si, and S2 
on £^ 

.51 



>f2 



Pi 



(65) 



/3' 



such that 

£Fif3 = I3'£f2j lF2<^ = 0i^Fi, Pa — id + silp^, a/3 = id + S2^f2- 

Example 7. Consider two presentations and ([M)) of the KdV equation from 
Example m The operators of diagram (l65t are: 



and 





— 6uDx 


— 6u 




-1 












A + 


6v 6u 





/3 = (1 0) 
si = 0, 



-A. -1 













S2 














1 






The form of operators a and /3 is obvious from the form of operators a and b in 
Example m The operators a' and /3' show how equations and ([M)) are obtained 
one from the other. 
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Example 8. Let the number of dependent variables m be equal to 1. Consider the 
lift L : J°° (tt) — !> J°° (vr) of the Legendre transformation 

It is defined wherever det ^ and preserves the Cartan distribution. Of 

course, L is not a map of fibre bundles. 

Consider an equation £ defined by F = 0, then L*{F) — defines the same 
equation £. The operators in diagram (I65p are as follows: 

Fi=F, 

F2^L*{F) 

and 

a = —1, a' = 1, si = 0, 
13 = -I, 13' = I, S2 = 0. 

Let us explain how to compute the maps a and /3. To this end, we have to take a 
symmetry of the Cartan distribution X G Xc{t^) /CX{tt) and see how the generating 
section of X transforms under the Legendre map. The generating section can be 
computed by the formula if = lu(X), where uj = du — w^-i dx^ is a Cartan form. 
So, a{ip) = uj{L{X)) = L*{uj){X) = -io(X) ^ -Lp. The same holds for /3. 

2.3. Symmetries and recursions. We have defined symmetries of a differential 
equation £ as elements of quotient space 

syuY£ ^ Xc{£)/CX{£), 

where Xc{£) = {X ^ X{£) \ [X,eX{£)] C CX{£) }, that is, symmetries of equa- 
tions are symmetries of the Cartan distribution on it modulo trivial symmetries 
(the ones belonging to the Cartan distribution). 

Obviously, symmetries of an equation form a Lie algebra with respect to the 
commutator. 

Given an inclusion £ C J°° (tt) , each symmetry X £ sym £ contains exactly one 
vector filed X^ S Xc{£) vertical with respect to the projection tToo : £ M, where 
M is the base manifold of the bundle tt. Next, for every such a field there 
always exists an evolutionary vector field E^> such that 

For any such a field, the restriction ip = ip'\^ £ >i = k >f(7r)|^ depends on X"" only. 
Hence, we can denote the field X" by E^. 

The element S >f such that X"" = E^ is called the generating section (or a 
characteristic) of this symmetry. 

If the equation at hand £ is defined by an equality F = 0, then the existence of 
a symmetry Eip boils down to the condition 

E^,{F%=Q, 

where, as above, ip' G ><(7r) is an arbitrary extension of (p, which is equivalent to 
the condition 

Ipi^p) =0 on (66) 
This is the determining equation for the symmetries of the equation £ = {F = 0}. 
The Jacobi bracket dTSl) yields a Lie algebra structure on symf in terms of 
generating sections: 

{(p,'ijj} = e^{ip) -i^iilj), (p,^; ekerip Oi- (67) 

Searching symmetries, that is solving (j66p . one usually begins with choosing 
internal coordinates on £. 
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Example 9. We have seen in Example S] that on the KdV equation 

Ut — 6uUx — Uxxx = 

the functions x, t, Us = Ux...x (s times) can be taken to be coordinates on £. These 
are internal coordinates for the KdV equation. 

Of course, the choice of internal coordinates is not unique. 

Next, one can start with finding symmetries whose generating sections depend 
on derivatives of order less than some number k, tp — ip{x, t,u,ui, . . . , Uk). 

Example 10. Let us find symmetries of the KdV equation such that Lp — Lp{x, t, u, ui). 
The determining equation for symmetries has the form: 

(A -Dl- QuDx - 6ux)ip{x, t, u, Ml) = 0. (68) 

The left-hand side of this equation is a polynomial in and U2 ■ The coefficient of 
the product M2U3 is equal to — 3(/9„j„j, so ipmui = 0, and we have 

ip — 'p^{x, t, u) + (p^{x, t, u)ui. 

Now, the left-hand side of (|68l) is a polynomial in U3, U2 and ui, with coefficient 
of M3 equal to —SDxif^)- Hence, 

ip = ip^{x, t, u) + (p^{t)ui. 

With such a (p, the left-hand side of does not depend on M3 any more. The 
coefficient of the product 1*1^2 is — 3<y9^„, hence 

if = p°^{x,t)+uip°\t,x)+ip'^{t)ui. 

The coefficient of U2 is ~3px^, so that 

(p = (^°° {x,t) + UlfP^ (t) + ip^ (t)ui 

and the left-hand side of (p5)) is a polynomial in ui and u. The coefficient of the 
product uui shows that tp'^^ = 0. The remaining coefficients of ui and u, and the 
free term reveal that 

where cq and ci are arbitrary constants. Therefore we have found two independent 
symmetries of the KdV equation 

fi = Ux and ip2 — QtUx + 1- (69) 

Let X e sym£ = Xc{^)/CX{E) be a symmetry. Vector fields Y G Xc{£) be- 
longing the equivalence class X we shall call representatives of X. As we explained 
above, any symmetry has a unique representative of the form E^. But this repre- 
sentative can be not the simplest one. 

Example 11. Symmetries ((69)) of the KdV equation can be represented, respec- 
tively, by the fields Yi = —d/dx and Y2 — d/du — 6td/dx, because 

Eu^-Yi = Dx, Eetu^+i-Y2 = 6tDx. 

The fields Yi and I2 are the lifts from the zero order jet space, hence they correspond 
to one-parameter groups of transformations; 

Yi: x' ^ X — e (translation along x) , 

Y2: x' = X — Get, u' — u + e (Galilean symmetry). 
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A symmetry X is called classical if it can be represented by a field Y\^, with 
Y £ X{tt) being the lift from a finite order jet space. Classical symmetries form a 
subalgebra of the Lie algebra sym£. By the Lie-Backlund theorem (see [l3]), Y is 
a lift from the zero order jet space (if the number of dependent variables ra > 1) or 
from the first order jet space (if the number of dependent variables m = 1). Thus, 
in coordinate language, the generating section of a classical symmetry has the form 




{ipi,. . . ,ipm) for m > 1, 
ipi^x"^ ,u,Ui) for m = 1, 



where (pa = ^0(2:% o-kix\u^)u'^, while ip is an arbitrary smooth function. 

The vector field Y on J° (tt) or (tt) that represents the symmetry with generating 
section ip has the form 

y^iET=ibj^-El=i<^k^ forTO>l, 

^v^n d^_d_,( spn dv_\ , sr^n (d^, d^-._d_ for rn - 1 

[ 2^i=l dui dx^ ^ l^i=l"'^du,) du^ l^i=iydx- ^ ^^du> dui ^Ol HI — 1. 

Example 12. In Example fTOl above we computed symmetries of the KdV equa- 
tion with generating sections depending on x, t, u, and u^. To find all classical 
symmetries, we allow also dependence on ut and compute symmetries in the same 
manner as in Example [TU] We get two additional classical symmetries: 

ifT, = Ut and ipi = xuj. + 3tut + 2u, 

with the corresponding one-parameter groups of transformations being 

933 : t' = t — e (translation along t), 

ip4: x' — e~^x, t' = e~^^t, u = e^^u (scale symmetry). 

Solving (j66p for sections that depend on at most fcth order derivatives will 
not give us a complete description of all symmetries. We can find all classical 
symmetries or a slightly larger subspace of symmetries, which can be considered to 
be a lower estimate of the full symmetry algebra. Letting the maximal order k be 
arbitrary and by solving (l66t describe the dependence of ip on derivatives of order 
fc. A; — 1, etc. This will be an upper estimate of the symmetry algebra. Sometimes 
these estimates allow to find more symmetries and, in some cases, obtain a complete 
description of the symmetry algebra. We refer the reader to [13] for examples of 
such calculations. 

A different approach is to look for a recursion operator that is a C-differential op- 
erator R: X ^ X such that there exists another C-differential operator R' satisfying 
the condition 

ipR = R'ip- (70) 

Operators of the form R = D^i?, with □ being arbitrary, enjoy ([70)) for all equations, 
so we consider recursion operators modulo such trivial ones. Obviously, i?(sym£) C 
symf , so that having a recursion operator we can produce an infinite number of 
symmetries from a given one. 

Example 13. The heat equation ut — u^x has two recursion operators of the first 
order 

R^ = and R2 = 2tD^ + x 
(and, of course, the identity operator, which is of no interest). 

For nonlinear equations we often are able to find a nontrivial recursion operator 
only if we allow it to contain the "integration" operator . In Section [331 we 
explain how to define such recursion operators in a rigorous way. 



26 



JOSEPH KRASIL'SHCHIK AND ALEXANDER VERBOVETSKY 



Example 14. As it was already mentioned, the KdV equation has the Lenard 
recursion operator 

R = Dl + 4u + 2ua;D-^. (71) 

It is obvious that R{ux) = Ut and R^{ux) = R{ut) = Uxxxxx + WuUxxx + '20uxUxx + 
30u^Ux- Thus, we have a new higher (that is non-classical) symmetry of KdV. We 
can proceed in the same way and compute R'^{ux), R'^{ux) and so on. As it was 
already mentioned (see Section [TT51) . all R'^lux) exist, that is the computation of 
will be possible for all k. So, we have constructed an infinite set of symmetries 
of the KdV equation. 

If we try to apply the Lenard recursion operator to the other two classical sym- 
metries, (^2 (Galilean) and (^4 (scale), we get R{'^2) — 2(^4, but Ri^fi) does not 
exist. In fact, R{(pi) is a nonlocal symmetry, as explained below in Section [31 

Now, let us explain how to compute recursion operators. To this end, note that 
C-differential operators k ^ >t modulo operators of the form \^lp can be naturally 
identified with elements of Aj(£) ^^^(f) >f, that is with valued Cartan 1-forms. 
This identification takes an operator R: >c -k to the form w/f G Ac(^^) ®T{£) ^1 
such that ujr{E^) = R{'^)- 

Next, recall that the Cartan 1-forms uj G Aj(f) arc functions on the tangent 
covering T{£) (see Example (H]) that are linear along the fibres of the projection 
r: T(f ) — J> f . Hence the forms lo G Aj(f) (8'jf(£) x are elements of the pullback k 
of the module x by r. Thus, to a recursion operator R: x ^ x we assign a 
fibre-wise linear element ujn G x. 



In coordinates, to the operator R 



X)/ /"'^-D/ll there corresponds the element 



with being coordinates along the fibres of r. 

Condition (j70p on R is equivalent to the following condition on un: 

ipiujR) =0. (72) 
T{£) 

Example 15. For the heat equation ut = u^x the tangent covering is given by 

Ut = Uxx, Vt = Vxx- 

Let t, Uk = Ux...xj and Vk = Vx...x be internal coordinates on it. To compute 
recursion operators of order 1 in Dx we are to solve the equation 

{Dt - Dl){f{x, t, Uk)vx + g{x, t, Uk)v) = 0. 

It is easy to show that the solutions are: v, Vx, and 2tvx+xv, with the corresponding 
operators id, Dx, and 2tDx + x, as it was indicated in Example 1131 

Example 16. The tangent covering over the KdV equation has the form 

Ut — Uxxx — GuUx = 0, Vt — Vxxx — GuVx — 6UxV = 0. 

Computation of recursion operators amounts to solving equations like the following: 

{Dt ~Dl~ 6uDx - 6ux){f2Vxx + fivx + fov + f-iV-i), (73) 

where fi — fi{x, t, Uk) and w_i is a new variable such that 

Dx{v-i)^v, (74) 

Dt{v^i) ^ Vxx + &UV. (75) 

As we saw above, the variables vi correspond to the total derivatives D/, so that 
the variable w_i is defined in ([7^ to correspond to D^^ , while ((75)) provides the 
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equality Dt{Dx{v-i)) — Dx{Dt{v^i)). We shall defer a rigorous explanation of this 
computation until Section [31 

One can check that v^x + 4uu + 2uxV^i is a solution of ()73p , it yields the Lenard 
recursion operator ((7T|) . 

The Lie algebra structure on "ksv Ip C x (the Jacobi bracket on symmetries) has 

a natural extension to a Lie superalgebra on ker £ p , called the Nijenhuis 

A2(f:)®^(£)^ 

(also Frdlicher- Nijenhuis) bracket and denoted by |- , •]. For a detailed definition 
we refer the reader to [S2J . Since we identified recursion operators with elements of 

ker ip , we can compute the Nijenhuis bracket of them. Two particular 

cases are of importance for us here: the bracket of a recursion operator with itself 
(the Nijenhuis torsion): 

^{R, i?,](<pi, <P2) = {Ri^i), R{iP2)} - R{R{vi), iP2} - R{ipi,Riip2)} + RH^u^2} 

= ■^fl,V2(-R(<Pl)) - ^R,vi{R{'P2)) + R^R.iptim) - RiR..V2i'P2), 

where ip2 G >c, and the bracket of a recursion operator and a symmetry with 
generating section (p £ keilp C x (the Lie derivative): 

L^{RW) = I^, RW) = R{^')} - RW, p'}, e X, 

or 

L^{R)^E^{R)^[£^,R]. 
A recursion operator R is called Nijenhuis (or hereditary) if its Nijenhuis torsion 
vanishes, i.e., i?] = 0. Almost all known recursion operators are Nijenhuis, 
including the operators encountered above. The main property of Nijenhuis opera- 
tors is the following: for every two symmetries with generating sections ipi, ip2 € x 
such that {ipi,ip2} — 0, L^-^{R) = 0, L^.-^(R) — and for arbitrary ki and k2 we 
have {R''^{(pi),R''^{^2)} = 0. 

Example 17. For the KdV equation take ipi = (p2 = Ux- Obviously, we have 
Lu^{R) — 0, where R is the Lenard recursion operator ((7T|) . so that all symmetries 
R''{ux) commute. 

2.4. Conservation laws. A conserved current u on an equation £ with n inde- 
pendent variables is a closed horizontal (n — l)-form on £, i.e., a form uj G Af^^^{S) 
such that dftW = 0. 

Example 18. Consider the equation of continuity in fluid dynamics 

Pt + [pv^)x^ + {pv^)x2 + {pv^)xs = 0. 
The form u! = p dxi A dx2 A dx^ — pv^ dt A dx2 A dxs + pv^ dt A dxi A dxs — pv^ dt A 
dxi A dx2 is a conserved current. 

Example 19. For the KdV equation ut — Uxxx — QuUx — the forms 

udx + {uxx + 3u^) dt, 

dx + (2uuxx -^1+ 4m^) dt, 

(u^/2 — u^) dx + [uxUxxx — u^j./2 — iu^Uxx + 6uu^ — 9u'*/2) dt 
are conserved currents. 

For any 77 G N^~'^{£) the form w = dur] is always a conserved current. Such 
currents are called trivial because they are not related to the equation properties 
and hence are of no interest. Consider the quotient space 

Hr\£) = {^^K'\£) I duLo = Q}/{^^kl~\8) \ ^^dHV. v&K-\£)} 
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of horizontal cohomology of £. We also want to quotient out the topological con- 
served currents that he in the image of the map C: H''-\£) H]^'^{£) induced by 
the natural projection A"~^(f ) — >■ AJJ~^(f ); here H^~^{£) is the (n — l)st group of 
the de Rham cohomology of the space £. Such currents are related to the topology 
of the equation £ only. Thus, we define c\{£) = HJ^^^{£)/ im( to be the set of 
conservation laws of equation £. 

We shall now discuss how to compute conservation laws for normal equations. 
To this end, let us consider the following complex: 

o^n%£) ^n\£) ^n^{£) ^ (76) 

where 

n'{£) = cm, ^"in = z^/ Mmn Ar'(f )), (7?) 

and ZP = kerd/i C A^{£) A^"^(f), for p > 0. The differential 5 is induced 

by the Cartan differential dc, so that 6^ = 0. 

Remark 17. Complex ([75)) is a part of Vinogradov's C-spectral sequence (see |1261 
IMlinniina), namely ri^^l^) - EP/'-\£) for p > 0, n°{£) = eI^''-\£) / H''-^ {£) , 
and 5 = d*'""^. 

Assume that f is a normal equation defined by a normal section F G P. In this 
case complex (j76p can be described in the following way: 

np{£) = ep/ep, (78) 

where Qp is a subset of C^i{>c, P) that consists of operators A e C^i{k, P) such 
that 

p-i 

e*pA{<fi,...,ipp^i) - ^ A*" {(fi, . . . , epi^pa), ■ ■ ■ , (fip-i) = 

for all iy9i, . . . , fp-i G ^, where *° denotes the operation of taking the adjoint with 
respect to the ath argument. The subset C Qp consists of operators A G 0^ of 
the form 

A(^i, . . . , ^j,_i) = ^(-l)"+iA'(£"j.(^„), ^1, . . . , . . . , ^p_i) (79) 

a=l 

for some C-differential operators A': Px>cx---x>c^P. 
In particular, for p = 1, we have 

n\£) ^keitp C P. (80) 

If p = 2 then 

n^i£) = {AeCi>c,P) I tpA = A*eF}/{A'£F \ A' eC{P,P), A'* = A'}. 

The differential S: n°{£) ^ n^{£) is given by the formula S{u}) = V*(l)|£, 
with V being a C-differential operator from P to A5J(7r) such that dhto — V(i^) 
on J°°(7r). 

The differential S: nP{£) nP+^{£), p>l, has the form 
p 

5(A)(^i, = 5](-l)"+'^A,^i,...,^„,...,^,(^a) + V|/^ (^1, . . . , ipp), 

a=l 

where V is a C-differential operator V : P x x[tt) x • • • x x{ti) iir(7r) that satisfies 
the relation 

p-i 

tpA{Lpi,...,Lpp_l) - ^ A*"- {ipi, . . . ,ip{Lpa), . . . ,ipp-l) = V{F,ipi,...,(pp_l) 

(81) 
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on J°°(7r). 

In the case of an evolution equation F = ut ~ f we have A G C^^*{>c, k) and the 
operator V can be chosen in the form 

V(/l,i^l,...,(Pp_l) = -•^A,¥;i,....Vp-i(^)' 

where h € P. 

If p = 1 and tp E imd C fl^{£) then we can put V = 

The above description of elements in flP{£) in terms of C-difFerential operators 
makes sense only for a given inclusion of the equation £ to a jet space. If we consider 
two inclusions 

J°°(^i) 



£ 




so that the corresponding linearizations are equivalent and we have diagram (|65l) . 
then the operators Ai and A2 that define the same element of VI'p{£) with respect 
to two inclusions are related as follows: 

Ai = a'*A2(a( ■),... ,a{-)), A2 = /3'*Ai(/3( ■),..., 

Elements of Q}[£) are said to be co symmetries of the equation £. We say that 
isomorphism (j80p takes a cosymmetry to its generating section (or the characteris- 
tic) belonging to ker£J, C P. 

The "Two- line Theorem" by Vinogradov (see, e.g., [13]) implies that complex (j76p 
is exact at the term ri°(f ), hence the conservation laws of £ form a subset of the 
space of cosymmetries, cl(£^) C Vl^{£). The generating section of a cosymmetry 
that belongs to 5{c\{£)) is called the generating section of the conservation law. 

As noted above, to compute the generating section of a conservation law we 
extend it arbitrarily to a form w G A'^~^{Tr) on the jet space J°°{tt), so that d^wj^ = 
0. Hence there exists a C-differential operator A: P ^ such that dhUi = 

A(F); the element V' = € ^ is the generating section of the conservation 

law under consideration. 

The generating section ^ of a conservation law can always be extended to the 
jet space J°°{tt) in such a way that {4>, F) = dhOJ, with w being a conserved current 
for the same conservation law. 

Remark 18. The above describe procedure is one of the ways to compute the differ- 
ential 5: r2°(f) f2^(f) in (ffS)) for an arbitrary equation £. If £ is presented 
in an evolutionary form then computation of the generating section is simpler 
and more straightforward. Namely, let t, x^, . . . ,x" be the independent variables 
and [uj] G ^l'^{£), where 

n 

w = X d.!^ A • • ■ A + ^ T; dt A dx^ A • • • A dV A • ■ • A dx". 
Then the corresponding generating section is 

where d/5u^ is the variational derivative with respect to . 

Example 20. The generating section of the conservation law from Example [T51 is 
equal to 1. 
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The generating sections of the conservation laws from Example [T51 are 1, 2u, 
and Uxx + 3m^. 

The determining equation for the conservation laws (or, to be more precise, of 
cosymmetries) of the equation £ — {F — 0} is 

i*FW = 0. (82) 

This equation is dual to equation ((M)) for symmetries. The computations involved in 
solving (|82p are very similar to those used to compute symmetries. However, unlike 
the case of symmetries, not all solutions of ((82)) give conservation laws. To check if 
the generating section of a cosymmetry corresponds to a conservation law we can 
use the following corollary of the "Two-line Theorem" by Vinogradov (see [I^): if 
the de Rham cohomology H'"-{£) — 0, complex (1751) is exact at the term n^{£). 
So, in this case the generating section ip e ker£p is the generating section of a 
conservation law if and only if S^ip) = 0, that is, there exists a self-adjoint operator 
A' = A'* e C{P, P) such that 

+ V|/ = A'ip, 

where V G C{m:{tt), P) satisfies the equality 

£*p{i;)^ViF) onJ°°(7r). 

Note that the generating section ?/> of a conservation law can always be extended 
to the jet space J°°(7r) in such a way that the horizontal n-form {ip,F) will be 
exact: {ip,F) = d/iCj, with uj\^ being a conserved current that corresponds to tp. 

2.5. A parallel with finite-dimensional differential geometry. III. In this 
section we begin compilation of a dictionary between the geometry of normal dif- 
ferential equations and finite-dimensional differential geometry, similar to the one 
for jet spaces from Sections [1 . 71 and [TT^ 

We start just as we did for jet spaces: we consider the space £ endowed with 
the Cartan distribution C and take for the points of that "manifold" the maximal 
integral submanifolds of C, i.e., the solutions of £. Further, the dictionary reads: 



Manifold M Normal equation £ 

points i — > solutions 

functions C°°(M) < — > conservation laws cl(f) 

the de Rham complex i — > complex ([75| 

of differential forms . . . ^ I^p-^ (f) ^nP{£)^ ... 

vector fields < — > symmetries 

the tangent bundle < — > the tangent covering r : T{£) £ 

In addition to considerations from Section 11.71 on jet spaces, this dictionary is 
justified by the following facts. 

First, on a finite-dimensional manifold the differential forms are functions on the 
tangent bundle with odd fibres. Correspondingly, definition (|77p shows that 

n*{£) = c\{T{£)). 

with elements of ^^{£) given by fibre- wise p-linear conserved currents. Fibres of 
the tangent covering r : T[£) — > £ are assumed to be odd. 

Since elements of r2^(f) can be understood as conservation laws on T{£), we 
can ask what are the generating section of these conservation laws? For the el- 
ement of ^'P{£) that corresponds to an operator A e C^_i{x,P) the generating 
section is (— V*%A), where the operator V is given by ([ST]) . Here we interpret 
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skew-symmetric C-differential operators >f x • • • x >f — ?> Q, modulo operators of the 
form ([75]) . as elements of the module Q pulled back on T{£). This describes the 
isomorphism (j78|) . 

Second, on a finite-dimensional manifold two natural actions of vector fields on 
differential forms exist: the interior product and the Lie derivative. Correspond- 
ingly, on an equation £ the evolution field E^, ip e symf , induces the interior 
product and the Lie derivative on n*{£) defined by (1771) : 

v: np{£) ^ np-^{£), L^: np{£) ^ np+\£). 

These operations are related to the differential 6 by the usual identity 

In terms of C-differential operators, the interior product : ilP{£) — > flP~^{£) 
for p > 1 is the contraction of the operator with ip. For p = 1, the interior product 
iipiip), Ip € kei £*p, arises from the Green formula and is the conservation law defined 
by the conserved current a;|^ e A^^^(£) such that 

V) - {v,epW) = dhOJ on J°°M. 

If G P is a generating section of a conservation law, then a symmetry (p £ sym £ 
acts on it by the formula 

with some operator □ G C(P, P) that satisfies the equality (-pif) — '^{F) on J°°{n). 

Third, on a finite-dimensional manifold a symplectic form gives rise to a Poisson 
bracket. The corresponding construction for an equation £ relies on the notion of 
a symplectic structure that is a closed element of r2^(f). We do not assume that 
the symplectic form is non-degenerate, so the Poisson bracket will be defined on a 
subset of c\£ (recall that conservation laws are analogues of functions on £.) 

In terms of C-differential operators, a symplectic structure is the equivalence 
class of operators A e C(>f, P) such that 

tpl\ = l\*iF, ^A,^.(¥'l)-^A,^i(<y52) + V|/H^1,(P2)=0, (83) 

where (pi, G V : P x >f — > >f is a C-differential operator such that 

- A*^F = V(P, • ) on J°°(7r), 
modulo operators of the form A7f, A' e C(F, F), A'* = A'. 
Example 21. For the simplest WDVV equation 

^yyy ~ ^xxy ~^ ^xxxUxyy = 

the operator is a symplectic structure. 

For evolution equations conditions (j83|) amount to 

A* = -A, eA,^Av2) - iA,^M = (-X^M)- 

The construction of the Poisson bracket is similar to the one on a finite-dimensional 
manifold. Let Q, € n^(f ) be a symplectic structure, i.e., S{n) = 0. A conservation 
law with the generating section tp is called admissible if there exists a symme- 
try ip G sym£ such that 

i^ = i^{n). (84) 
Symmetries that correspond to admissible conservation laws in the sense of ([M)) 

are called Hamiltonian symmetries. 

By definition, the Poisson bracket of two admissible conservation laws lj and w' 

with the generating sections '0 and "0', respectively, has the generating section 

{w,w'}o = Lipi'^'') = iipiip'^, 
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where Lp and ip' are Hamiltonian symmetries corresponding to conservation laws 
w and w'. To the conservation law {wjw'jsi there corresponds the Jacobi bracket 
lyj'}, so that Hamiltonian symmetries form a Lie algebra: [L,^, L^pi] = L^^ ,^/y. 
As we see from (|84p. a symplectic structure takes a symmetry to a conservation 
law, in terms of operators (j83p . this map has the form — A{lp). 

2.6. Cotangent covering to a normal equation. In the previous section we 
discussed geometry related to the tangent covering and functions on it (forms). 
But what about the cotangent covering? 

We have defined the tangent covering for an equation £ without fixing an in- 
clusion of f to a jet space. Dualizing such an invariant definition requires use of 
rather complicated homological algebra, so we will define the cotangent covering 
for an equation £ embedded to a jet space, £ C J°°(7r), and then check that the 
construction does not depend on the choice of the embedding. 

For a normal equation £ given by a normal section F = 0, F £ with P being 
a module over J°°{n), we define the equation T {£) C J^{P) by the equalities 

F = 0, iUp) = 0, 

where the tilde denotes the pullback to (P) and p £ P corresponds to the iden- 
tity operator P P under the identification P = C{P, P). The natural projection 
T* : T {£) ^ £ is called the cotangent covering to £. 

In coordinates, we havep — (p^, . . . if the coordinates on J^{P) are x*, Uj,pj, 
with and being fibre coordinates along projections £ ^ M and J^{P) £, 
respectively. 

Below we assume that not only equation £ is normal but the operator £p is 
normal (see Remark as well. Then T* {£) will also be a normal equation. 

Remark 19. Obviously, for every £ the cotangent equation T*{£) is an Euler- 
Lagrange equation with the Lagrangian density L = {F,p). In applications, con- 
sidering T {£) instead of £ is occasionally useful to handle the equation as though 
it were Lagrangian (see, e.g., [351 Volume 1, Sections 3.2, 3.3]). We refer to [371 
Section 4.5.1] and |105l Section 5] for more details and references. 

If we have two inclusions of equation £ to jet spaces 

J°°(^i) 



£ 




then the adjoint linearizations t*p_^ and fp^ are equivalent: 



Pi 

A 

* 

Y 



A 



(85) 



}i2 



such that 



l\a'* = a*i%, tpX = f3*tp^, a'*P'* = id + sl'tp^, p'*a'* = id + s^l^, 
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where the operators a, a', j3' , si, and S2 are defined in (|55|). Therefore, the 
cotangent coverings constructed using i*p_^ and l*p^ are isomorphic, thus the cotan- 
gent coverings do not depend on the choice of inclusion £ C J°°(7r). 

Now we describe an isomorphism between symf and the subspace of clT {£) 
that consists of conservation laws with fibre-wise linear conserved currents. Thus 
we will justify the first parallel in the prolongation of our dictionary: 

Manifold M Normal equation £ 

the cotangent bundle i — > the cotangent covering t* : T* {£) — )> £ 
multivector fields i — > conservation laws cl(T*(f)) (86) 

Let <p & xhe the generating section of a symmetry of £. Extend it to an element 
Lp S xr(7r) and consider the Green formula 

{If{v). V-) - tp{^)) = V), (87) 

where tp ^ P, uj{ip,'4') G A'^^^^{Tr). The mapping ijj i— > uj{ip,ip) is a C-differential 
operator P — > AJ^^"^(7r), so that it gives rise to a closed form € A'^~^{T* {£))■ 
The induced map sym£ — >■ clT'*{£), which takes the symmetry with the gener- 
ating functions (f to the conservation law with the current w^, gives the desired 
isomorphism. 

In fact, formula (1571) gives more. It holds not only for generating sections of 
symmetries of £, but also for an arbitrary tp E >c, so that we obtain a map x — ^ 
A'^~^{T* {£))■ Since >ir is a direct summand in >c{T*{£)), the element uj{ip, tp) yields 
an element p G ^1^{T {£))■ One can prove that p does not depend on the choice of 
inclusion £ — > J°°(7r) used in its construction. 

In terms of isomorphism ([75)1 . we have p = {p, 0). 

This is a very important element since it plays the role of the canonical 1-form 
pdg on a finite-dimensional cotangent space: 

Normal equation £ 

pen\r{£)) 

canonical symplectic structure 

n = 5{p) e n^T* {£)) 

In terms of operators (j78p . the canonical symplectic structure on T* {£) has the 
form 

-(° 0.) 

As to entry (I86p of our dictionary, we take it for the definition of multivectors. 
Since we are interested in skew-symmetric multivectors, we assume the fibres of the 
cotangent covering T* {£) £ iohe odd. 

We call conservation laws of T {£) whose currents are fibre-wise p-linear vari- 
ational p-vectors on £ and denote their set by Dp{£). Thus, Dq{£) = cl£ and 
Di{£) = sym£. 

For Dp{£) we have a description in terms of C-differential operators similar 
to (HH]). Namely, 



Manifold M 

the canonical 1-form p dq 

the canonical symplectic form 
dp A dq 
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where is a subset of C^i{P, k) that consists of operators A G C^i{P, k) such 
that 

p-i 

£fA(V'i, • • • , V-p-i) - E (^1' • ■ • ' • ■ • ' '/'p-i) = 

a=l 

for aU ■01, . . . , V'p-i G P, where *° denotes, as before, the operation of taking adjoint 
with respect to the ath argument. The subset C Sp consists of operators A S Sp 
of the form 

ACV'l, . . . , Vp-l) - ^(-l)"+lA'(£^(l^„), Vl, . . . , V-a, . • • , V'p-l) 

Q = l 

for some C-differential operators A' : k x P x ■ ■ ■ x P ^ x. 
In particular, for p — 2, we have 

D2{£) = { A e C(P, ^) I If A = A*tp }/{ A'tp \ A' e C{k, x), A'* = A' }. 

For the element of Dp{£) that corresponds to an operator A G C^i{P, x) the 
generating section is (— V*^, A), where the operator V:PxPx---xP^-P is 
given by 

p-i 

ipA^il^i, ^p^i) - ^ A*° (V-i, . . . , ^^(V'a), . . . , V'p-i) = V(F, Vi, . . . , V-p-i)- 

(88) 

In the case of evolution equation F = ut — f we have A G C^^*{k, k) and the 
operator V can be chosen in the form: 

V(/i, Vi, ■ • • j-iAp-i) = ^A>i,...>p_i(M> 

where h € P. 

Since we assume the fibres of cotangent covering to be odd, the bracket defined 
on cl(7~ {£)) by the canonical symplectic structure will be the variational Schouten 
bracket 

I,-] - DkxDi^Dk+i-i. 
In terms of C-differential operators, this bracket has the form: 

[Ai,A2l(V'i,...>fe+;-2) = (-l)"^A„v.,i.-i,(^i(V'<T(/,fc+i-2))) 

^ (-l)-A2(Vr(V.(i,.)),V'.(fe+i,fc+;-2)) 

+ (-1)'"' E (-irAi(vr(^.(i,o),V'.(m.feH-i-2)), 

where Ai G Dk{P), A2 G Di{P), Vi and V2 are defined by ([HHl) , V'l , • ■ ■ , ^/'/c+i-2 e 
P, cf. gH). 

The above description of variational multivectors in terms of C-differential op- 
erators makes sense only for a given inclusion of equation f to a jet space. If we 
consider two inclusions 
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SO that the corresponding hnearizations and adjoint Hnearizations are equivalent 
and we have diagrams (|55|) and ([55)) . then the operators Ai and A2 that define the 
same element of Dp{£) with respect to the two inclusions are related as follows: 

A2 = aAi(a'*( •),...,«'*(•)), Ai = /3A2(/?'*( • ), . . . ,/3'*( • )). 

An element A G D2{£) is called a Hamiltonian structure on £ if [A, A] = 0. 
Two Hamiltonian operators Ai and A2 are said to be compatible if |Ai,yl2] = 
(cf. Section [TT 



Remark 20. Using this definition of a Hamiltonian structure, the Hamiltonian for- 
malisms on jet spaces, including the Magri scheme, explained in Section II. 8[ can 
be extended straightforwardly to the case of equations described here. 

Example 22. The Camassa-Holm equation [T7] 

has a bi-Hamiltonian structure: 

Ai = Dx, A2 = -Dt - uDx + Ux- 
This equation is often written in the form 

mt + umx + 2uxm = 0, 
m — u + Uxx — 0. 
Then the bi-Hamiltonian structure takes the form 



A'^i °] A' ^ ° 

1 ^ Dx - Dl r ^ \ 2mDx + nix 



Example 23. Let £ he a bi-Hamiltonian equation given by F = and Ai and A2 
be the corresponding Hamiltonian operators. The Kupershmidt deformation [761 
ES] £ oi £ has the form 

where w = {w^ , . . . ,w^) are new dependent variables. For example, the KdV6 
equation [ST] 

Vt + Vxxx + l'2vVx -Wx=0, Wxxx + ^VWx + ^WVx = 0, 

is a Kupershmidt deformation of the KdV equation corresponding to the Hamil- 
tonian operators Ai = Dx and A2 = D^. + 8vDx + 4:Vx- 

The following two variational bivectors define a bi-Hamiltonian structures on £: 

iF+Al(w)+A*(w)} ' (y~^F+Al(w)+A*{w) 

Example 24. The equation 

Zyy + i^/z)xx + 2 = 

associated with an integrable class of Weingarten surfaces [7] is bi-Hamiltonian with 
operators D^ and 2zDxy — ZyDx + ZxDy. 
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3. Nonlocal theory 

Nonlocal phenomena in the theory of integrable systems are quite common. Here 
by nonlocality we mean an extension of the initial system by new variables (fields) 
that are related to the old ones by differential relations. Perhaps, the simplest 
way to observe how nonlocal objects originate is to analyse the action of recursion 
operators on symmetries. 

Example 25. Consider recursion operator ((55)) R = D^ + 4:U + 2uiD^^ that gener- 
ates the higher KdV equations (it can be shown that successive application of R to 
the first symmetry ipi = ui results in polynomial expressions in m, mi, . . . , w^, . . . , 
see, e.g., [571 lll3j V When one applies the operator R to the first (x, i)-dependent 
symmetry 

1 

ifl =tui + - 
6 

(the Galilean boost), this will result in the scaling symmetry 

1 N 2 
<P3 = tU3 + [6tu + -x)ui + -u, 

but application of the recursion operator to (^3 leads to an expression that contains 
the nonlocal term D^^{u) which can not be expressed in the geometrical terms 
introduced aboveQ. 

An apparent way to incorporate this nonlocal object into the initial geometric 
setting is to introduce a new variable, say w, that is related with the old one by Wx = 
u. This relation, due to the KdV equation, implies another one: wt = + u^x 
and thus we shall result in the system 

Ut = 6uUx + Uxxx, Wx=U, Wt^Su'^ + Uxx- 

A general geometric formulation of this construction was first introduced in [1331 
171] and below we shall give a concise exposition of the theory together with a number 
of applications. 

3.1. DifTerential coverings. The notion of a covering was already used in Sec- 
tion [2] in the context of the tangent and cotangent coverings. Here we discuss it in 
more detail. 

Let 6 C J°°(7r), where tt: E ^ AI, dimM = n, be an equation. Consider a 
locally trivial bundle r : f — > f and endow the manifold £ with an n-dimensional 
distribution C in such a way that 

(1) C is integrable and 

(2) for any point 6 £ £ the restriction drj^^ is a one-to-one correspondence 

between Cg and the Cartan plane C^^gy 

Then we say that ip is endowed with the structure of a differential covering (or 
simply a covering^ to be short) over £. 

Coordinates. Consider a trivialization of the bundle r and let , . . . , . . . be 
fibre- wise coordinates (the so-called nonlocal variables). The number r of these 
coordinates is called the dimension of r. 

Let Di, . . . , Dn be the total derivatives on £. By Property (0) in the definition 
of the distribution C there exist r- vertical vector fields Xi , . . . , X„ on £ such that 
the fields 

bi = D.I+ Xi, i = 1, . . . , n, 



Of course, the Lenard operator itself contains a nonlocal summand, but we can consider it 
just as a convenient reformulation of the Magri relation I I54II . 
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lie in C. Then Property ([T|) is equivalent to the system of equations 

D,{Xj) - Dj{Xi) + [X,,Xj] =0, l<i<j<n, (89) 

where Xi,...,X„ are r- vertical fields and Di{Xj) denotes the component-wise 
action. Since the vector fields Xi are r-vertical fields, they can be presented in the 
form 

where Xf are smooth functions on £, while £, as a manifold with distribution, is 
isomorphic to the infinite prolongation of the system of PDEs 

— z = j = l,...,r, 

which extends the initial equation £ and is compatible over it due to This 
system is called the covering equation. 

Example 26. Consider the one-dimensional covering over the KdV equation de- 
termined by 

d d 

ow ow 
The covering equation in this case is 

dw dw 2 

ox at 

and is isomorphic to the potential KdV equation wt = ?>w^ + Wxxx- 

Note the the relation between w and u may be expressed in the form w ~ J u dx, 
or w — D~^u and thus this is exactly the nonlocality that arose in Example 1251 

Example 27. Let again the base equation be the KdV and the covering be de- 
scribed by the system 

X = u + w^+X, T ^U2 + 2wui+2u^ + 2{w^ - \)u-4:\{w^ + 1), (90) 

where A G M. Actually, (j90p determines a one-parameter family of covering struc- 
tures in the trivial bundle f x M — > but the covering equation is isomorphic to 
the modified KdV equation wt = Qw^Wx + Wxxx for any A. Of course, the covering 
under consideration is a geometric realization of the Miura transformation |,95^ . 

Remark 21. Note that any differential substitution u — ip{x,'w, . . . ,wi, . . .) is 
associated with a covering over the initial equation, though this covering may 
be infinite-dimensional. For example, such is the covering over the KdV equa- 
tion itf — 6uUx + Uxxx = determined with the Hirota substitution 

u = -2—- Inw 
OX'' 

(see [5D]). Nevertheless, in spite of the infinite dimension of this covering, the 
covering space is isomorphic to the fourth-order scalar equation 

WWxt - WtWx + WxxxxW - ■iWxxxWx + 3w^^ = 

in one unknown function. 

Example 28. Let P be the module of sections for some vector bundle ^ over £. 
Then the bundle jj^ : Jf^{P) £ oi horizontal jets is an infinite-dimensional cov- 
ering over £. If are adapted coordinates in J^{P) then the total derivatives 
lifted to J^iP) are of the form 

(91) 
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This construction is generahsed in the next example. 

Example 29 (A-coverings). Let E be an equation and consider a C-differential 
operator A: P Q, where P and Q are modules of sections for some vector 
bundles ^ and C over Z. Let $a : Jh°{P) ^ Jh°{Q) be the corresponding mor- 
phism of vector bundles (see Proposition II. 3|) . Then, under natural conditions of 
non-degeneracy, £a = ker $a is a sub-bundle in ^^o '■ J/^ (P) £ that carries a 
natural structure of a covering: the total derivatives in this covering are obtained 
by restriction of the operators (l9T|) to £a- We call this covering the A- covering 
over £. 

If the operator A is locally given in the matrix form A = || J^k '^afs^KW then 
the subspace £a C J/^iO is described by the relations 

-iK 



= 



.K 



and their prolongations. Obviously, the tangent and cotangent coverings are par- 
ticular cases of this construction with A — ig and A = ig, respectivelja. 

A-coverings play the key role in solving the following factorisation problem: 
let A' : P' ^ Q' be another C-differential operator; how to find all operators A: P ^ 
P' such that 

A'ov4 = PoA, (92) 

i.e., such that the diagram 

P Q 



P' ~^Q' 

is commutative? Note that any operator A of the form A = B' o A, where B' : Q ^ 
P' is an arbitrary C-differential operator, is a solution to Such solutions will 

be called trivial. 

To find nontrivial solutions, first note that since A' is a C-differential operator it 
can be lifted to the covering just by changing the total derivatives Di to the lifted 
ones Di. Denote this lift by A'. Second, let us put into correspondence to any 
operator ^ = || J2k '^apW ^^'^ vector-function 



.K a,K 



, r' — dimP', 



Then one has: 

Proposition 3.1. Classes of solutions of Equation (192p modulo trivial ones are in 
one-to-one correspondence with solutions of the equation 

A'i^A) = 0. 

Operators satisfying (|92p take elements of ker A to those of ker A'. 



Consider system (j89|) that determines a covering structure in the space £ and 
assume that the coefficients of the vertical vector fields Xi are independent of 
nonlocal variables w". In this case, (j89p reduces to 

D,{Xj)^Dj{X,), l<i<j<n; (93) 



Here and below we consider normal equations and use the notation instead of £p which is 
justified by the results of Subsection 12.21 
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the corresponding covering is called Abelian. The covering in Example [52] is an 
Abelian one, while the covering associated with the Miura transformation (Exam- 
ple [27l) is not. 

Let dim t = 1 and define a differential horizontal f-form on £ by setting 

n 

UJr = '^Xidx\ (94) 
i=l 

Then ([55]) amounts to the equation 

dhior = 0, (95) 

where d^ is the horizontal differential on £. Thus, one-dimensional Abelian cover- 
ings over £ are in one-to-one correspondence with closed horizontal {n — l)-forms. 

In Example 1271 we presented a one-parameter family of coverings over the KdV 
equation. Are these coverings different for different values of the parameter A or not 
and what the word "different" means in this context? The answer is the following. 

Consider an equation £ and two coverings Ti: £i ^ £ , i — 1, 2, over £. We say 
that these coverings are gauge eguivalent (or simply equivalent) if there exists an 
isomorphism ip: £i ^ £2 of the equations £1 and £2 such that the diagram 



£1 




£2 



£ 



is commutative, i.e., T2 o ip = n. In this sense, all coverings (j90p are different, 
i.e., pair- wise non-equivalent for different values of A. A general cohomological 
technique to check whether a parameter can be eliminated or not was suggested 

in [iniina- 

We say that a covering t: £ £ is trivial if for any point 9 & £ there exists a 
neighbourhood U 3 9 such that 

(f ) is a trivial bundle; 

(2) there exists an adapted coordinate system in U for which the fields Di are 
of the form Di — Di, i = 1, . . . , dim Af . 

Theorem 3.2. There exists a one-to-one correspondence between equivalence classes 
of one-dimensional Abelian coverings over £ and elements of the horizontal coho- 
mology group Hj^ {£) given by (|94p . In particular, a covering t is trivial if and only 
if the form ujr is a co-boundary, i.e., LOr — d/i(/). 

Since Hl{£) coincides with the group of conservation laws when dimAf — 2, 
Theorem 13.21 allows one to construct special type of coverings by conservation laws 
of the equation at hand. 

Example 30. The Camassa-Holm equation 

ZUxtlxx '^'^xxx 

admits the conservation law 

cj = (u - Uxx) d.T + ^(m^ - iv? + 2uuxx) dt; 
consequently, the corresponding covering is given by 

Wx Uxx, wt = ^{ul - -I- 2uUxx)- 
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Remark 22 . Since the group Hj^ {£) is trivial for normal equations in the case dim M > 
2, this implies that such equations possess no nontrivial one-dimensional Abelian 
covering. Actually, there exist very strong indications that these equations do not 
have finite-dimensional coverings at all (see [87]). 

Example 31 (the KP equation). Consider the dispersionless Kadomtsev-Petviashvili 
equation 

It admits an obvious covering 

Wx^Uy, Wy ^ Ut - 6uUx + Uxxx, (96) 

which at first glance seems to be one-dimensional. But this is not the case, because 
Equations (I96p do not contain information on the derivative wt- To incorporate 
these data, we must introduce infinite number of nonlocal variables w^, w^,... 
such that 

1 r r+1 

W — W, Wf. — W , . . . , Wf. — W , . . . 

and express their x- and y-derivatives using (j96p . Thus, the covering is infinite- 
dimensional actually. 

It was shown above that one-dimensional Abelian coverings can be constructed 
using conservation laws of the equation. Another type of coverings is related to 
Wahlquist-Estabrook prolongation structures |135[ 11361 [55] and their description is 
based on the following ansatz: Let ut = f{u,ui, . . . ,Uk) be a system of evolution 
equations, u = (m^, . . . , m™), / — (/^, . . . , /™) being vectors and Ui denoting the 
ith derivative with respect to x. Let us look for coverings such that the coefficients 
of the fields X and T in 

Dx^Dx+X, Dt = Dt+T 

depend on u, mi, . . . , Uk-i and nonlocal variables only. Then, locally, the description 
of such coverings locally reduces to representations of a certain free Lie algebra (the 
so-called Wahlquist-Estabrook algebra) in vector fields on the fibre W of the trivial 
bundle t: £ xW ^ £. 

Example 32. Consider the potential KdV equation 

Ut ^ ul+ Uxxx (97) 

and let us describe coverings Dx ^ Dx + X , Dt = Dt + T over £ such that the 
fields X and T depend on u, ui and U2 only. Straightforward computations show 
that all these coverings are of the form 

X = ti^a + ub + c, 

T = (2uu2 - uj + 2u'^ui)a+ {u2 + 2uMi)b + ui [c, b] + iu^[b,d] +u[c, d] + e, 

^ (98) 

where a, b, c, d and e are vector fields on the fibre W of the covering (i.e., such 
that they do not depend on the equation coordinates) which enjoy the commutator 
relations 

2a=[a, b], b=[a, c], d = 2c+[c,b], 
[a, d] = [c, e] = 0, 

[b, d] + ^ [b, [b, d]] = 0, [b, e] + [c, [c, d]] = 0, 
[a,e] + [b,[c,d]] + i[c, [b,d]] =0. 
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Now, to find all Wahlquist-Estabrook coverings for (|97p amounts to describing 
representations, as vector fields on W, of the free Lie algebra generated by the 
elements a, b, c, d, and e with the above relations. 

If W — M then all such representations, up to an isomorphism, are 

d d d 

aK^— , h>~^{2w + l3) — , {w^ + I3w + -f) — , 

ow ow ow 

ow ow 
where /3, 7 G M and A = — 47. The corresponding one-dimensional coverings, 
up to gauge equivalence, are of the form 

X = (u'^ + 2wu + w'^ + 7)t;— 

ow 

(the parameter /3 can be removed by a gauge transformation) and with T given 
by (|98p : they are pair- wise inequivalent for different values of 7. 

Remark 23. The term covering also refers to the parallel between classical dif- 
ferential geometry and geometry of PDEs. Namely, if we define dimension of an 
equation £ (or of a jet space J°°(7r)) as that of the corresponding Cartan distribu- 
tion (i.e., the number of independent variables) then fibres of a differential covering 
become zero-dimensional, and this complies with the definition of a topological 



covering. That was the initial reason to name the object in |133| . 

But the parallel goes far beyond this trivial observation. In [52] , a new power- 
ful invariant (the fundamental Lie algebra) of differential equations was proposed 
whose role in the theory of differential coverings is quite similar to the one that the 
fundamental group plays in topology. In particular, the fundamental Lie algebra al- 
lows one to enumerate (locally) all coverings over a given equation in the same way 
as conjugacy classes of subgroups of the fundamental group enumerate topological 
coverings. So, the dictionary evolved in the previous sections can be continued: 

Manifold M Differential equation £ 

topological dimension i — differential dimension 
topological coverings i — differential coverings 
fundamental group < — > fundamental Lie algebra 

Note also that using the fundamental Lie algebra technique the author of [52] 
proved nonexistence of Backlund transformations for some pairs of differential equa- 
tions. It seems that it is impossible to achieve such a result by other methods. 

3.2. Nonlocal symmetries. The concept of a symmetry discussed in Section [5] 
can be generalised to the nonlocal situation. Consider an example. 

Example 33. Let 

Ut = UUx + Uxx (99) 

be the Burgers equation (its Lie algebra of symmetries was fully described in |133j ). 
Direct computations show that ([M]) does not possess symmetries of the form ip = 
if{x, t, u), but if one extends the setting by a new (nonlocal) variable w such that 

Wx = u, Wt ^ + Ux (100) 
then the equation ig (ip) = will acquire a new family of solutions of the form 

ip^ {au'^2ax)e-^'", (101) 
where a = a{x,t) is an arbitrary solution of the heat equation at = Uxx- 
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The question is: can functions (jlOip be considered as symmetries of Equa- 
tion (|M1) in some natural sense? To answer this question, consider an arbitrary 
equation £ C J°°(7r) and a covering r: £ — > f. We say that ip is a, nonlocal sym- 
metry (or r-symmetry) of £ if it is a symmetry of £. 

Coordinates. Let £ C J°°{-k) be an equation and t: £ — ^ £ be a covering locally 
given by the total derivatives 

+ z = l,...,dimA^. 

j 

Then any nonlocal r-symmetry is of the form 

3 

Here ip is an m-component vector-function on £ that satisfies the equation 

hiv) = 0, (103) 

are functions on £ such that 

b.,{^^) = ij^) + Y.^r (104) 



and 



a 

over internal I 
coordinates 

(recall that the "tilde" over a C-differential operator denotes its natural lifting to 
the covering). 

Example 34. Let us consider Example [33] again. In the case of covering (llOOp 
Equations (|104p take the form 

D^ii:)^^, DtW=u^ + D,{^). (106) 

Consequently, for ip of the form (|10ip we see that 

= -2ae-5» 

satisfies (|106l) . Thus, the pair of functions ip and iJj determine a nonlocal symmetry 
of the Burgers equation in the sense of the above definition. 

However, the situation of the previous example is not generic. 

Example 35. Consider the covering 

Wx = U, Wt = -|- Uxx 

over the KdV equation ut = Quu^ -\- Uxxx and let us try to find nonlocal symmetries 
in this covering. In the case under consideration, Equations (jl04p acquire the form 

Dx{'ip)^V, Dt{ij)^6uv + Dliip), (107) 

while pII3)) is 

Dtip) = Guiip + 6uDx{p) + Dliip). 
The simplest solution of the last equation that depends on w is 

p — tu^ + ( lOtu + — a; ) U3 + 4 ibtui + — ) U2 + 2 ilbtu^ + xu + -rw ) ui + —u^- 



But solving (I107P with p) of the above form leads to a contradiction: no function -0 
exists on £ such that (|107p is valid for our p. 
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Nevertheless, if we introduce another nonlocal variable w' satisfying 



u 



then (|107p will be resolved in the new setting. 

But a similar problem arises at the next step: now we need to reconstruct the 
coefficient ip' at d/dw'. 

The procedure we encountered in Example [35] is typical and we shall describe 
it in general terms now. Let r : £ — ^ £ be a covering and denote by T and T the 
function algebras on £ and £, respectively. An M-linear map X : 7^ — > is called a 
T-shadow if 

(1) X is a derivation, i.e., 

X{fg) = fXig) + gXif) 

for all /, 5 e J"; 

(2) the action of X preserves the Cartan distribution, i.e., Lx{oj) G Ac(^) as 
soon as w G Ac{£) (or, equivalently, for any Cartan field Y on £ and its 
restriction Y to the commutator [X, Y] XY — YX is a Cartan field 
again) . 

In particular, any symmetry of £ can be considered as a shadow in an arbitrary 
covering r. 



Coordinates. Let U be the set of internal coordinates on £. Then any r-shadow is 
given by the formula 

ouj 



u]eu 

where (p^ , . . . , ip"^ are functions on £ and Di , . . . , £)„ are total derivatives on £ 
(cf. (jlOSp ). The vector function ip — {ip^ , . . . , 1^9™) must satisfy the equation 

££{ip)=0. 

We say that a r-shadow X is reconstructed in r if there exists a nonlocal r- 

symmetry X such that X = X. As Examples [Ml and [55l show, not all shadows 

can be reconstructed in a straightforward way. A general result that describes the 
reconstruction procedure was proved in [BB] (see also [7T]): 

Proposition 3.3. Let t: £ ^ £ be a covering and X be a T-shadow. Then there 
exists another covering f: £ ^ £ and a f -shadow X such that X\j- — X . 

Thus, putting t — tq and t,;+i = fi and applying Proposition [331 sufficiently 
(maybe, infinitely) many times we shall arrive to a covering in which the given 
shadow is reconstructed. 

Coordinates. Actually, the results of [BB] not just state the existence of the needed 
covering but provide a canonical way to construct the one. The construction is in 
a sense tautological and mimics relations (|104[) . 

Proposition 3.4. Let t: £ ~^ £ be a covering over an equation £ with nonlocal 
coordinates , . . . , , . . . and Di , . . . , -D„ be total derivatives in this covering. Let 
also ip be a r-shadow. Then: 

(1) the relations 

^=4H^)+Ea^*"' * = 1,---,", J = l,...,dimr, (108) 

a 

deRne a covering over S whose dimension equals that of r; 
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(2) equations (|104p are solvable in this covering. 

Remark 24. It follows from (|108p that for an Abelian covering r the covering f 
is Abelian as well. Hence, at every step of reconstruction obstructions to solv- 
ing (|104p lie in the horizontal cohomology group of the corresponding equation. 
Consequently, ii t: £ — >■ f is an Abelian covering and Hl{£) — then any r- 
shadow can be reconstructed to a nonlocal r-symmetry. In particular, any local 
symmetry of £ can be lifted to £. 

Let £ be an equation and {[w"]}, G A\{£), be an R-basis of the group Hl{£). 
Assume that 

uj" = dx^ + ■■■+ dx" 
and consider the covering ti: £i ^ £ determined by 

for all a and i = l,...,n. For £i let us construct the covering T2- £2 — >■ £1 in 
a similar way, etc. The covering t^: £,, £ obtained as the inverse limit of the 
sequence 



is called the universal Abelian covering over £. By construction, Hj^{£^,) = 0. 

Proposition 3.5. For an arbitrary Abelian covering t: £ £, there exists a 
uniquely (up to a gauge equivalence) defined morphism 

c f 

£* ^£ 



£ 

and any r-shadow can be reconstructed in r* . In particular, any symmetry of £ 
can be lifted to f * • 

Remark 25. Though the covering f whose existence is stated in Proposition 13.31 is 
determined canonically by the shadow X, the new shadow X is not unique, but is 
defined up to an infinitesimal gauge symmetries of t, i.e., up to F S sym{£) such 
that Y\jr — 0. Due to (|104p . these symmetries are given by the equations 

DM')=Y.d^.r (109) 

a 

and are of the form Y = ''P" d / dw" . 

Example 36. For Example E51 Equations (jl09l) take the form 

L),(V') = 0, AW=0 

and consequently infinitesimal gauge symmetries are ^d/dw in this case, 7 S R. 
On the other hand, if we consider Example [57] then Equations (|109p are written 

as 

D^{ij)^2wip, Dtiijj) = 2{ui + 2uu)-4:Xvu)tJj. 
The only solution of this system is "0 = and thus there is no ambiguity in shadow 
reconstruction in this case. 

Non-uniqueness of the solution to the problem of reconstruction leads, in turn, 
to the problem of commutation for shadows: no well defined way to compute the 
Lie bracket of shadows is known. This problem was first indicated in |102j . A way 
to solve it was suggested in |125) , but a practical realization of the approach is 
somewhat cumbersome. 
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To conclude this subsection, let us make a remark also related to the problem of 
reconstruction. Let r : £ —i' £ he a. finite-dimensional covering and X be a symmetry 
of £. One can (at least, locally) lift X to f in an arbitrary way. Then, if X is an 
integrable vector field (i.e., if it possesses the corresponding one-parameter group 
of transformations) then the lifted field X is integrable as well. If X is a symmetry 
of £ then this means that we managed to reconstruct X up to a nonlocal symmetry 
in the covering r. 

Conversely, consider the one-parameter group of transformations {Ax} corre- 
sponding to X and for any A G M define an n-dimensional distribution on £ 

by 

C^:e.^C,^-iA,*(c^-i(,)) (110) 

where 9 £ , Cg is the Cartan plane at the point and denotes the differential 
of the map F. 

Proposition 3.6. Correspondence (jllOp determines a one-parameter family t\ of 
pair-wise inequivalent coverings over £ such that tq = t. 

Example 37. Take the covering 

X — u -h w'^ , T — U2 + 2wui 2v? + 2uP'u 

over the KdV equation and apply ProDOsition l3. 41 using the Galilean boost tui -1-1/6 
for the symmetry X. This will result in the Miura covering described in Example [?7l 

Not all one-parameter families of coverings can be obtained by this procedure (a 
counter-example can be found in f22j[21]). But a weaker result was proved in [54] : 

Theorem 3.7. Let t\: £ ^ £ be a one-parameter family of coverings regarded 
as a deformation of the covering t = tq. Then the corresponding infinitesimal 
deformation is a r-shadow. 

3.3. Backlund transformations and zero- curvature representations. In con- 
clusion, let us briefly discuss how the constructions of Backlund transformations |109l 
I108j and zero-curvature representations [T] are translated to the geometrical lan- 
guage. 

A Backlund transformation between two equations £' and £" with the unknown 
functions u' and u", respectively, is another equation £ in unknown functions 
both u' and u" such that for any solution u' of £' a solution u" of £ is a solu- 
tion of £" as well and vice versa. If £' coincides with £" then one speaks about 
auto- Backlund transformation. 

Example 38. Consider the sine-Gordon equation 

Uxy — sinu. (Ill) 

Then the system 

v — t— u 2 v — u 

- u,y = 2A sin — - — , -\- = - sin — - — , (112) 

Z A 2 

where A 7^ is a real parameter, determines the classical auto-Backlund transfor- 
mation (a one-parameter family, actually) for (jllip . see p7] . 

Example 39. The second example (which now can also be considered as a classical 
one) was found in [135j . It is of the form 




(113) 
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and relates solutions of the KdV equation to each other (or, to be more precise, 
system (jll3p is a Backlund transformation for the potential KdV equation, while 
solutions of the KdV itself are obtained hy u — Vx)- 

Analysis of these two examples (as well as other ones) shows that a Backlund 
transformation between equations £' and £" is a diagram 



£ 




£' £\ 

where r' and r" are coverings. The correspondence between solutions of £' and £" is 
achieved in the following way. Let v! — u'{x') be a solution of £' and assume that r' 
is a finite-dimensional covering. Then the Cartan distribution of the equation £ 
induces on the finite-dimensional manifold = (t')~^(u') C £" an n-dimensional 
integrable distribution. In the vicinity of a generic point the latter possesses a 
(dimT')-parameter family of maximal integral manifolds that are projected to u' 
by t' and to the corresponding family of solutions of £" by r". Generically, such a 
correspondence is non-trivial provided r' and r" are not gauge equivalent. 

Example 40. A common way to construct non-trivial Backlund transformations 
is the following. Let t: £ ^ £ he a. covering and / : f — > f be a finite symmetry 
of f , i.e., a diffeomorphism preserving the Cartan distribution. Then the com- 
position r' = r o / is a covering as well and the pair (t, t') is an auto-Backlund 
transformation for £. If / is not a gauge equivalence then this transformation is 
non-trivial. 

Consider covering (|90)) from Example [27] and note that the change of the nonlocal 
variable w O —w is a symmetry of the covering equation (the mKdV one), but is 
not a gauge symmetry of the covering itself. Thus, for any value of the parameter A 
we get an auto-Backlund transformation of the KdV equation, i.e., a one-parameter 
family of Backlund transformations. 

Note that the Wahlquist-Estabrook construction (Example [51]) is a consequence 
of the latter one. 

Remark 26. Families of Backlund transformations, like the ones from Examples 1381 
and 1391 give one an opportunity to construct special exact solutions of integrable 
equations (such as multi-kink solutions for the sine-Gordon equation, multi-soliton 
solutions for the KdV, etc.). The construction uses the nonlinear superposition 
principle which, in turn, is based on the following informal statement: 

Theorem 3.8 (the Bianchi Permutability Theorem). Assume that an equation £ 
possesses a one-parameter family of auto-BUcklund transformations B\ and let A 6 
R be the parameter. For any solution u — u{x) of £ denote hy B\{u) the set of 
solutions obtained from u hy means of B\. Then for any Ai ^ A2 there exists a 
solution uai,A2 G B\-^{B\2{u)) n Bx^lBx^^lu)) that is expressed as a bi-differential 
operator appUed to some sohitions ui G B\-^{u) and U2 G B\^{u). 

This "theorem" was first observed by Bianchi in [TU] (see also |109j ) in applica- 
tion to the sine-Gordon equation (Example [55)1 and since then dozens of examples 
were computed, but nevertheless a general formulation of this statement (and, con- 
sequently, its general proof) is unknown to us. Some hints to a rigorous approach 
to the problem can be found in [89j . 

Geometrical theory of Backlund transformations is also related to an unortho- 
dox approach to recursion operators [55]. Consider an equation £ and its tangent 
covering t: T{£) £ (see Section [2]). Recall that symmetries of £ arc identified 
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with sections of r that take the Cartan distribution on 8 to that on T{£). Hence, 
if we consider a diagram of the form 



T{£) T{£) 




where r' and r" are coverings, then this Backlund transformation wih relate sym- 
metries of £ to each other. Thus, this Backlund transformation plays the role of a 
recursion operator for symmetries of £. 

Example 41. Consider the KdV equation ut = QuUx + Uxxx and two copies of its 
tangent covering with the new dependent variable v that enjoys the the additional 
equation 

Vt = QUxV + QuVx + Vxxx- 

Introduce a nonlocal variable v by setting 

Vx = W, = Qui) + Vxx- 

Thus, internal coordinates in £ are 

X, t, U^Uq, Ux = Ui, . . . , V = Vo, Vx = Vi, . . . , V. 

Define the covering r' by 

r': {x,t,Uk,Vk,v) ^ {x,t,Uk,Vk) 

and the covering t" by 

r': {x,t, Uk, Vk,v) H> {x,t,Uk, D^{v2 + 4uw + 2uiv)). 

The Backlund transformation obtained in such a way is the geometrical realization 
of the Lenard recursion operator R ~ + Au + 2uiD~^ given in Equation ((55|) . 

Another, less trivial example will be considered later (see Example H51 below) . 
after discussing the concept of zero- curvature representations (ZCR). 
Let T : 5 — > 5 be a covering. We say that it is linear if 

(1) r is a vector bundle; 

(2) the action of vector fields Di, . . . , Dn on J-{£) preserves the subspace of 
fibre-wise linear functions. 

Coordinates. Let v^, . . . , w'', ... be local coordinates along the fibre of r and the 
covering be given by the total derivatives 

A=A + 5]xr^, z = l,...,n. (114) 

r 

Then the covering is hnear if and only if the coefficients in fjll4p are of the form 

x: = E^>"' 

a 

where X^^ are smooth fimctions on £. If we now identify the vertical terms Xi = 
Xld/dv^ in (jll4p with the function- valued matrices 

X,. = 



48 



JOSEPH KRASIL'SHCHIK AND ALEXANDER VERBOVETSKY 



then (|114p will be rewritten as 

bi ^ D.I+ Xi, i = 1, . . . , n, 

while the conditions [Di, Dj\ =0 will acquire the form 

D,{Xj) - Dj{X,) + [X„ X,] =0, l<i<j<n. 

In other words, we arrive to the classical definition of a ZCR (cf. [T]). 

Example 42. The well known two-dimensional ZCR for the KdV equation (see [I]) 
is given by 

D,^D^ + A, bt^Dt+ B, 

where 

^ ^ / 1\ ^^f -u^ 2ii - 4A 

\-u + X Oy ' - 2?i^ + 2Xu + 

Example 43 (vacuum Einstein equations). Consider the Lewis metric ds^ = 
2/(x, y) dx dy + X]j<j 9ij i'^ with coordinates x, y, z^, and (see [7S]). 

The the vacuum Einstein equations read 



(i/det gg^g ^)y + i\/detggyg ^ 0- (115) 
After a re-parameterisation, (|115p acquires the form 

(116) 



- lUx +Uy _ VxU + y + U^Vy IVx+Vy 



u 2 x + y ' " u 2 x + y 

Backlund transformations and ZCR for (jll6p were constructed in many papers (see 
e.g., [illlllMllHS]). The latter is of the form 



D,=D.,+ A, Dy = Dy+ B, 



where 



^^2[ie-l)v., i^j' ^=2^1(1-^)., l^j' 



here 9 — •\/(A + y)(A — a;) and A is the spectral parameter. 

Using ZCR (|117l) . a three-dimensional covering over T{£) can be constructed 
(see [HI]). Let [/ and V be the variables in T{£) corresponding to u and w, respec- 
tively, and , up' , w'^ be the nonlocal variables. Then the covering is given by the 
relations 

1 l-e 2 1 + 61 3 1 + 6*^^ 1 + 9 

= —z—v^w + ——v^w — Ux + —pu.^lJ, 

2 2u'^ 2u 2u'^ 

2 1 + 9 1 1 + 61 2 1 + 6* rr l+^T. 
W^= 5— WxW u^w ^'"xU + -—T-Vx, 

wl^{9- l)v^w^ + ——Uxw'^ -I —Vx 



and 



1 61-1 2 1+^' 3 1 + ^ TT 1 + 



3 l-6» 1 1 + 6* 3 1-61^^ 
- -^-y^ + -J^^y^ + -^Vy 

This covering gives rise to a Backlund transformation of the form 

9U' = 2uw^ + U, 9V' = -u'^w^ -w^, 

i.e., to a recursion operator for symmetries. 



1 + 61 1 1+61 2 1+6* rr 1 



v 
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Remark 27. Note that with an arbitrary covering r; £ — > £ one can naturally 
associate a linear covering : £ ^ £. The space T" 5 is a submanifold in T £ 
and consists of tangent vectors that vanish under the action of the differential r^,. 
Note also the existence of the exact sequence of coverings 

^ T" £ ^ r£ ^ T" £ ^ 



where r'* : 7"** f — ^ f is the quotient. The terms of this sequence possess the 
following characteristic property: integrable section^Q of t" are infinitesimal gauge 
symmetries of r, integrable sections of f are, as it was mentioned above, nonlocal 
T-symmetries of £, and integrable sections of r" are r-shadows. 

Concluding remarks 

We described a geometrical approach to partial differential equations which 
proved to be efhcient both from the theoretical point of view and in a lot of appli- 
cations. Based on this approach, in particular, Hamiltonian formalism for arbitrary 
normal (2-line) equations is constructed. On the other hand, a number of interest- 
ing and important problems are waiting for their solution. We intend to continue 
the research along the following lines: 

• Generalisation of the Hamiltonian formalism from normal equations to ar- 
bitrary p-line ones that, in particular, include gauge-invariant systems. 

• Incorporation of Dirac structures into the above described scheme and elab- 
oration of their computation and use. 

• Further development of the nonlocal theory and, in particular, analysis of 
differential coverings over the systems with the more than two independent 
variables and generalisation of the theory of variational brackets to nonlocal 
structures. 
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